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Abstract 

This  study  used  simulation  and  experiment  to  design  and  test  foil  packets  for  use 
in  the  Activation  Foil  Integrated  Detection  System  (AFIDS).  The  initial  plan  to  activate 
foil  packets  outside  with  the  pulse  reactor  at  White  Sands  Missile  Range  (WSMR)  was 
not  possible  due  to  WSMR  not  having  safety  approval  to  take  the  reactor  outside.  As  an 
alternative,  the  concept  of  using  liquid  nitrous  oxide  inside  a  reactor  to  simulate  large 
volumes  of  air  was  investigated. 

Simulation  using  the  Standardized  Computer  Analyses  for  Licensing  Evaluation 
(SCALE)  program  was  used  to  select  ten  foils  to  be  included  in  the  foil  packet.  The  size 
was  selected  with  a  target  activity  of  25  Bq  for  each  foil  four  days  after  the  activation 
when  located  500  meters  from  a  five  kiloton  equivalent  nuclear  weapon. 

We  analyzed  whether  N2O  could  replicate  large  volumes  of  air  in  neutron 
transport  experiments  since  one  cubic  centimeter  of  liquid  N?0  contains  as  many 
molecules  as  1371  cubic  centimeters  of  air.  A  neutron  propagating  through  the  liquid 
N2O  should  react  like  a  neutron  propagating  through  1371  times  as  much  air.  Actual 
testing  could  not  be  completed  at  the  Ohio  State  University  Research  Reactor  due  to 
hazardous  levels  of  radiation  in  the  facility  during  testing.  The  results  of  the  simulation 
showed  that  pragmatic  experimental  factors  inhibit  use  of  the  reactor  in  this  way  and  an 
outside  test  is  still  required  to  provide  the  information  Oak  Ridge  National  Laboratory 
(ORNL)  needs  to  validate  the  methodology  used  in  Domestic  Nuclear  Event  Attribution 
(DNEA). 
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MANUFACTURE  AND  TESTING  OF  AN  ACTIVATION  FOIL  PACKAGE  FOR 


USE  IN  AFIDS 

I.  Introduction 

Both  parties  in  the  2004  Presidential  Campaign  Debates  indicated  that  the  most 
dangerous  threat  to  the  United  States  was  a  domestic  nuclear  event.  Although  no  one 
likes  to  dwell  on  this  scenario,  we  must  be  prepared  to  deal  with  one.  Information  about 
the  device  involved  has  not  only  political  implications  but  is  important  in  relief 
operations  and  public  safety. 

Background 

A  nuclear  event  requires  analysis  to  detennine  weapon  type  and  characteristics  for 
use  in  Domestic  Nuclear  Event  Attribution  (DNEA).  Much  of  this  information  can  be 
gained  from  the  prompt  neutron  spectrum  of  the  weapon.  Materials  near  the  detonation 
will  be  activated  by  the  prompt  neutrons  and  infonnation  about  the  spectrum  can  be 
gained  from  the  analysis  of  these  activated  materials. 

One  difficulty  with  using  activated  materials  from  near  the  detonation  site  is  that 
the  exact  composition  of  the  materials  is  unknown.  The  proposed  Activation  Foil 
Integrated  Detection  System  (AFIDS)  is  a  passive  radiation  detection  system  consisting 
of  foil  packages  containing  foils  of  known  composition.  The  foils,  located  near  the 
detonation,  become  activated  by  prompt  neutrons  from  the  nuclear  detonation.  Since  the 
foils  would  be  chosen  to  activate  with  neutrons  at  certain  energy  thresholds  and  the 
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reaction  cross-sections  are  well  known,  the  activity  of  the  foils  collected  after  a  nuclear 
detonation  would  provide  fluence  and  energy  information  about  the  neutron  spectrum. 

The  neutron  flux  at  a  certain  point  can  be  calculated  by  computer  adjoint 
calculations  using  a  known  neutron  flux  from  a  known  source.  In  the  case  of  an 
unknown  source,  if  enough  information  is  available  from  threshold  activation 
measurements,  it  should  be  possible  to  unfold  the  neutron  source  spectrum  using  adjoint 
calculations.  This  method  has  been  proven  in  and  out  of  the  laboratory  in  limited  cases. 
Oak  Ridge  National  Laboratory  (ORNL)  has  developed  a  means  for  comparing  an 
unknown  spectrum  to  known  spectra  by  comparing  calculated  activities  from  the  known 
spectra  using  adjoint  fluxes  to  measured  activities  of  strategically  placed  foil  packages. 
After  a  nuclear  event,  the  activated  foil  packages  would  be  collected,  the  neutron  induced 
activity  measured,  and  the  results  sent  to  ORNL.  ORNL  would  then  compare  the 
measured  activity  to  activity  adjointly  calculated  from  the  known  weapon  spectra.  The 
comparison  would  rank-order  the  closest  matches  to  known  spectra  and  eliminate 
weapons  and  types  that  do  not  match.  The  features  of  the  most  closely  matched  spectra 
would  indicate  types  (e.g.  uranium  vice  plutonium)  and  help  to  validate  results  from  other 
methods  such  as  fractionation  analysis  [22], 

Problem  Statement 

The  purpose  of  this  thesis  is  to  design,  manufacture  and  test  activation  foil 
packages  by  simulation  and  testing  in  order  to  facilitate  further  research  that  provides 
several  measured  activities  for  ORNL  to  complete  proof-of-principle  verification  for  the 
methodology  to  be  used  in  DNEA. 
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Scope 


The  original  scope  of  this  thesis  was  to  design  and  build  foil  packages,  irradiate 
them  using  the  pulse  reactor  at  White  Sands  Missile  Range  (WSMR),  and  compare  the 
known  spectrum  found  by  adjoint  calculations  at  ORNL.  The  WSMR  pulse  reactor  does 
not  have  the  safety  approval  for  outside  testing  so  the  scope  was  modified  to  investigate 
the  concept  of  using  liquid  nitrous  oxide  (N20)  to  replicate  large  volumes  of  air  in  lieu  of 
outside  testing  at  WSMR.  Foil  packages  were  designed,  built  and  activated  and  neutron 
transport  through  N20  was  compared  to  neutron  transport  through  air  by  simulation  and 
experimentation  at  the  Ohio  State  University  (OSU)  research  reactor. 

Assumptions/Limitations 

The  foils  used  are  subject  to  physical  and  practical  constraints.  The  atomic 
content  of  the  foils  has  to  be  known  so  ORNL  can  do  pre-calculations  and  time  is  not 
spent  after  retrieval  determining  the  composition  of  the  foils.  The  reactions  of  interest 
have  to  cover  all  parts  of  the  neutron  spectrum  to  include  thennal,  epithennal  and  fast 
neutrons.  The  cross-section  of  the  reaction  of  interest  has  to  be  well  documented  to  allow 
ORNL  to  perform  calculations.  To  prevent  self-absorption,  the  mean  free  path  of  the 
resulting  gamma  must  be  greater  than  the  material  thickness  through  which  the  gammas 
will  be  measured.  For  the  purposes  of  packaging,  the  choice  of  foils  was  restricted  to 
ones  that  are  available  in  wire  form. 

Reactions  of  interest  from  neutron  activation  were  limited  to  ones  that  produce 
gamma  rays  during  decay  to  facilitate  counting.  Detectors  used  to  measure  the  gamma 
activity  of  the  foils  should  have  resolution  better  than  10%  full  width  at  half  maximum 
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(FWHM)  and  1%  absolute  efficiency  at  1  MeV.  For  design  purposes,  1%  absolute 
efficiency  was  assumed  and  0.1  counts  per  second  as  the  lower  limit  of  counting  based  on 
counting  with  a  High  Purity  Germanium  (HPGe)  detector  for  two  hours.  Foils  were 
assumed  located  500  meters  from  a  frve-kiloton  nuclear  weapon  detonation  and  the 
activity  measured  four  days  after  the  nuclear  event  for  design  purposes.  Since  the  lower 
limit  of  counting  is  0.1  counts  per  second  and  assuming  100%  decay  percentage,  100% 
branching  ratio,  and  1%  absolute  efficiency,  the  minimum  activity  of  the  foil  four  days 
after  activation  must  be  at  least  10  counts  per  second.  The  maximum  activity  of  the  foil 
four  days  after  activation  was  limited  to  1000  counts  per  second  to  minimize  pile-up  and 
dead  time  and  to  allow  the  same  geometry  for  all  foil  counting.  Higher  activities  can  be 
counted  by  placing  the  foil  farther  from  the  detector  but  the  system  must  then  be 
recalibrated. 

To  locate  a  peak  in  gamma  spectroscopy,  the  software  used  in  this  thesis  uses  the 
second  differential  to  locate  a  peak,  fits  a  gaussian  to  the  peak,  and  provides  the 
difference  between  data  and  background  counts  in  standard  deviations  [6],  A  standard 
deviation  of  three  was  assumed  as  the  minimum  detectable  peak. 

The  foil  packages  are  also  subject  to  physical  and  practical  constraints.  They 
were  designed  to  survive  a  five-kiloton  nuclear  detonation  at  500  meters.  This  requires 
surviving  8  psi  overpressure,  925°C,  shock  and  blast.  Since  the  exact  location  of  the 
blast  is  unknown,  the  foil  packages  cannot  be  oriented  towards  the  blast  and  must  be 
designed  to  capture  fast  and  epithermal  neutrons  from  any  direction.  Direction  is  not  as 
important  for  thermal  neutrons  because  thermal  neutrons  reach  equilibrium  in  velocity 
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and  will  backscatter  into  the  foils.  The  foil  packages  should  not  present  a  hazard  to  the 
general  population  or  the  environment.  Self-shielding  within  the  packaging  of  the  foils 
needs  to  be  accounted  for  and  minimized  so  each  foil  gets  maximum  exposure  to  the 
prompt  neutrons.  The  cost  per  foil  package  needs  to  be  low  since  a  significant  number  of 
foil  packages  are  needed  to  cover  potential  targets  and  low  cost  reduces  pilfering.  A  limit 
of  five  dollars  was  assumed. 

For  the  purposes  of  design,  the  simulated  point  source  was  monoenergetic  at 
14  MeV  (upper  energy  limit  of  bomb  neutrons)  and  yielded  10“  neutrons,  approximately 
the  equivalent  of  a  five-kiloton  weapon  [14].  The  duration  of  neutron  production  from  a 
nuclear  device  is  on  the  order  of  microseconds  while  the  half-lives  of  the  decay  elements 
are  on  the  order  of  hours  or  longer.  Thus,  the  decay  of  the  radionuclides  during 
activation  was  assumed  negligible  and  ignored  when  calculating  activation. 

Organization 

A  more  in-depth  analysis  of  the  problem,  design,  simulation  and  testing 
procedures  is  discussed  in  the  following  chapters.  Chapter  II  describes  the  theory  behind 
nuclear  weapon  neutron  output,  foil  activation  and  the  methodology  used  by  ORNL  to 
compare  neutron  spectra.  Chapter  III  provides  the  foil  package  design  procedures  and 
results  while  Chapter  IV  contains  the  testing  procedures  and  results.  Chapter  V 
investigates  the  concept  of  using  N?0  to  replicate  large  volumes  of  air.  Chapter  VI 
contains  conclusions  and  recommendations  for  future  study.  The  appendices  provide 
supporting  calculations,  tabulated  data,  and  references  used  throughout  the  thesis. 
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II.  Theory 


Definitions 

The  following  definitions  are  presented  in  the  context  that  they  are  used  in  this 

thesis. 

Foils. 

Foils  are  the  individual  wires  in  the  foil  packet.  They  are  chosen  because  of  their 
ability  to  capture  the  neutrons  in  the  energy  range  of  interest. 

Foil  Packet. 

Foil  packets,  or  foil  packages,  are  sets  of  foils  packaged  together.  Foil  packets  are 
the  physical  entity  deployed  in  the  field. 

Detectors. 

For  the  purposes  of  this  thesis,  detectors  refer  to  the  devices  used  to  measure  the 
gamma  rays  being  emitted  from  the  activated  foils.  Detectors  used  in  this  thesis  must 
have  a  minimum  absolute  efficiency  of  1%  and  minimum  resolution  of  10%  full  width  at 
half  maximum.  Detectors  are  not  to  be  confused  with  foil  packets. 

SCALE. 

Standardized  Computer  Analyses  for  Licensing  Evaluation  (SCALE)  software  is 
developed  and  maintained  by  ORNL  under  contract  with  the  U.S.  Nuclear  Regulatory 
Commission  and  Department  of  Energy  to  perform  criticality  safety,  radiation  shielding, 
spent  fuel  characterization,  and  heat  transfer  analyses  of  nuclear  fuel  facilities  and 
transport/storage  package  designs  [23],  For  the  simulations  in  this  thesis,  the  Shielding 
Analysis  Sequence  No.  4  (SAS4)  was  implemented.  SAS4  is  a  control  module  that 
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allows  calculation  of  radiation  doses  exterior  to  a  transport/storage  cask  via  a  three- 
dimensional  Monte  Carlo  analysis  [27].  The  output  of  SCALE  is  the  product  of  neutron 
flux  and  a  response  function.  Normally  the  response  function  is  the  factor  to  convert 
neutron  flux  to  radiation  dose.  In  this  thesis,  the  response  function  was  set  to  one  when 
the  required  output  was  the  neutron  flux  and  set  to  the  cross-section  of  the  reaction  of 
interest  when  the  output  was  used  for  calculating  activation.  The  fractional  standard 
deviations  in  the  output  of  SCALE  are  good  indicators  of  whether  the  simulation  is  good 
or  not.  The  goal  is  to  get  under  10%  fractional  standard  deviation  [4].  This  can  usually 
be  achieved  by  adjusting  the  number  of  particles  in  the  batches  or  the  number  of  batches 
used  in  the  Monte  Carlo  analysis. 

Effects  of  Nuclear  Weapons 

A  nuclear  detonation  is  similar  to  a  conventional  explosion  in  that  a  large  amount 
of  energy  is  released  within  a  limited  space  and  time.  This  energy  release  drives  the 
temperature  and  pressure  up  converting  all  materials  present  to  hot,  pressurized  gasses. 
These  gasses  expand  quickly  and  produce  a  shock  or  blast  wave.  A  nuclear  detonation 
differs  from  a  conventional  explosion  in  that  much  more  energy  is  released  in  a  shorter 
time  with  less  mass.  Much  higher  temperatures  are  achieved  and  more  energy  is  emitted 
as  heat  and  light  (thermal  radiation).  For  example,  a  foil  package  500  meters  from  a 
ground  burst  of  a  5  kT  device  will  experience  an  overpressure  of  8  psi  and  50  cal/cnL  of 
thennal  radiation  [8].  The  50  cal/cnr  of  thermal  radiation  translates  to  an  increase  of 
900°C  in  silicon  carbide  assuming  the  radiation  is  absorbed  in  the  first  millimeter  of 
material.  A  nuclear  detonation  also  emits  radiation  called  initial  nuclear  radiation  or 
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prompt  radiation.  The  radioactive  substances  left  after  the  nuclear  detonation  continue  to 
emit  radiation  that  is  known  as  residual  radiation  [8]. 

Nuclear  fission  (breaking  apart  of  atoms)  and  fusion  (combining  of  atoms)  both 
release  high-energy  neutrons.  A  pulse  of  neutrons  with  energies  of  approximately 
14  MeV  and  lower  is  emitted  from  the  detonation  and  is  called  the  prompt  spectrum. 
These  neutrons  activate  the  foil  packages  based  upon  their  energy. 

Radiation  diffuses  as  it  propagates  outward  from  a  point  source.  The  direct 
(unscattered)  fluence  is  given  by 


0  = 


N 

And2  ’ 


(1) 


where  N  is  the  number  of  neutrons  and  d  is  the  distance  from  the  source  [12].  Thus,  a 
point  source  of  1024  neutrons  (e.g.  a  5  kT  device)  would  have  an  unscattered  fluence  of 

14  2  13  2 

approximately  8x10  n  cm'"  at  100  meters,  3x10  n  cm'"  at  500  meters  and 
1.3  x  1012  n  cm'2  at  2500  meters. 


Activation 

Activation  foils  work  using  the  physics  of  neutron  absorption  in  isotopes  to 
produce  radionuclides.  These  daughter  radionuclides  then  decay  by  emitting  alpha  or 
beta  particles,  or  gamma  rays.  The  emitted  particles  or  rays  can  then  be  detected  and  the 
activity  of  the  foil  can  be  detennined.  The  activity  can  be  calculated  from 

A=m0,  (2) 

where  X  is  the  decay  constant  and  N0  is  the  number  of  activated  atoms. 
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The  change  in  the  number  of  activated  nuclei  is  the  rate  of  production  minus  the 
rate  of  decay.  Mathematically  this  is 

dN 

— — —  =  (fxrN T  —  AN ,  (3) 

dt 

where  </>oNT  is  the  rate  of  production  and  AN  is  the  rate  of  decay  or  activity.  The 
number  of  radioactive  nuclei  is  represented  by  N,  (f>  is  the  flux  of  neutrons,  cr  is  the 
absorption  cross-section,  Nr  is  the  number  of  target  nuclei  and  A  is  the  decay  constant 
[28].  The  absorption  cross-section  is  material  and  energy  dependent.  A  thermal  material 
has  a  large  absorption  cross-section  in  the  thermal  region  and  a  threshold  material  has  a 
large  cross-section  or  resonance  at  specific  energies  above  the  thermal  region  (0.5  eV). 
Therefore,  a  thennal  reaction  of  interest  would  be  initiated  by  thermal  neutrons  and  a 
threshold  reaction  of  interest  would  require  a  neutron  with  an  energy  corresponding  to  its 
threshold  to  initiate  the  reaction. 

The  pulse  length  of  neutrons  from  a  bomb  is  much  less  than  a  second  and  the  half- 
lives  of  the  radionuclides  of  interest  are  in  the  range  of  hours  to  years  so  the  decay  during 
activation  is  negligible  and  can  be  ignored.  Thus,  the  number  of  radionuclides  produced 
from  the  pulse  can  be  calculated  by 

N0=noNr,  (4) 

where  N0  is  the  number  of  radionuclides  produced  and  n  is  the  number  of  neutrons 

incident  on  the  target.  Because  no  activation  occurs  during  decay,  the  decay  can  be 
accounted  for  by 

N(t)  =  N0e-Xt,  (5) 
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where  t  is  the  time  after  activation  [13]. 

When  a  material  with  an  absorption  cross-section  is  exposed  to  a  constant  neutron 
flux,  the  induced  activity  builds  with  time  and  approaches  a  saturation  activity  Ajnf  as 
shown  in  Figure  1 .  In  this  example,  neutron  irradiation  has  proceeded  for  a  time  tG  when 
the  foil  is  removed  from  the  neutron  flux  with  an  activity  A0.  This  leads  to  a  saturation 
activity  given  by 


An  f  = 


(1 -«“*) 


(6) 


Once  removed  from  the  neutron  flux  the  radioactivity  will  decay  by  Equation  5. 


Figure  1 .  Activity  of  an  activation  foil  where  t0  is  the  time  of  the  foil's  removal 

from  the  neutron  flux  [12]. 

In  the  case  of  a  nuclear  detonation,  a  variety  of  materials  with  absorption  cross- 
sections  will  activate  with  exposure  to  the  neutron  flux.  The  composition  of  these 
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materials  is  unknown  so  the  starting  point  for  activation  is  unknown.  The  elements  used 
in  the  AFIDS  foils  may  be  common  in  many  materials  but  might  not  be  present  in 
combinations  or  concentrations,  may  be  shielded  so  as  to  mask  parts  of  the  spectrum,  or 
would  require  separation  techniques  which  reduce  or  eliminate  the  accuracy  of  the 
measurements. 

Neutron  Detection 

Radiation  detection  typically  requires  conversion  of  radiation  energy  into 
electrical  energy.  However,  neutron  detection  differs  from  other  detection  methods 
because  the  neutrons  have  no  charge.  Therefore  thennal  neutron  detectors  rely  upon 
neutron  induced  nuclear  reactions  to  produce  a  secondary  charged  particle  that  can  then 
be  counted  by  charged  particle  detectors  [30].  Although  these  detectors  can  detect 
neutrons  by  producing  a  prompt  output  pulse,  they  are  not  practical  for  AFIDS  as  the  cost 
of  the  detector  alone  prohibits  mass  fielding.  They  are  also  manpower  intensive  so  even 
if  the  detector  survived  the  effects  of  a  nuclear  detonation,  the  operator  might  not. 

Neutron  measurements  can  also  take  advantage  of  the  radioactivity  induced  in 
some  materials  by  neutron  interactions.  Foils  made  of  materials  selected  for  their  ability 
to  activate  with  neutrons  can  be  exposed  to  a  flux  of  neutrons,  collected  later,  and  then 
the  induced  radioactivity  can  be  counted  using  conventional  means  such  as  gamma 
spectroscopy.  This  passive  method  meets  AFIDS’  needs  of  low  cost,  low  maintenance, 
nuclear  event  survivable,  and  delayed  counting. 
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Gamma  Detection 


Two  common  types  of  detectors  found  in  the  lab  for  measuring  gamma  rays  are 
inorganic  scintillators  (e.g.  Nal(Tl))  and  semiconductor  detectors  (e.g.  HPGe).  The  two 
main  factors  in  choosing  between  them  are  energy  resolution  and  counting  efficiency. 
Energy  resolution  is  a  measure  of  the  ability  to  discern  fine  detail  in  the  incident  energy 
of  the  radiation.  Absolute  counting  efficiency  is  the  ratio  of  gamma  rays  detected  to  the 
total  gamma  rays  emitted  by  the  source.  It  is  dependent  on  the  geometry  of  the  setup  and 
the  energy  of  the  gamma  rays  as  well  as  the  detector.  In  comparing  the  two  detectors, 
Nal(Tl)  detectors  tend  to  have  better  counting  efficiency  while  HPGe  detectors  have 
better  resolution  [12].  Normal  absolute  efficiencies  for  HPGe  detectors  of  800  keV 
gamma  rays  are  around  one  percent  [7,  19].  Absolute  efficiencies  of  Nal(Tl)  detectors 
are  typically  about  an  order  of  magnitude  better  [12].  Typical  resolution  of  an  HPGe 
detector  is  1.95  keV  FWHM  at  1.33  MeV  and  the  normal  resolution  of  an  Nal(Tl) 
detector  is  about  an  order  of  magnitude  worse.  Better  efficiency  gives  more  counts  while 
the  better  resolution  may  let  multiple  foils  be  counted  simultaneously.  Either  type  of 
detector  should  work  to  measure  AFIDS  foil  packages  and  should  be  available 
throughout  the  country  at  colleges  and  universities  if  provision  and  planning  are  made  for 
their  calibration  and  use. 

The  detector  used  in  this  research  was  a  Canberra  Model  GC 10021.  The 
Canberra  Model  GC  10021  HPGe  detector  is  a  p-type,  closed-end,  coaxial  detector  that 
detects  gamma  activity  from  low  activity  foils.  The  crystal  has  dimensions  of  83 
millimeters  diameter  and  84.5  millimeters  length.  The  detector  was  configured  so  that 
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the  front  crystal  face  was  located  5  millimeters  from  the  detection  window  [30].  The 
detector  collects  electrical  charges  created  at  either  boundary  of  the  semiconductor 
material,  and  the  spectroscopy  system  records  these  pulses  and  differentiates  between 
energies.  The  detector  was  located  in  a  model  747  Canberra  Lead  Shield.  The  shield 
blocks  high  background  gamma  radiation  by  using  a  four-inch  thick  lead  wall  and  a 
1.0-millimeter  tin  and  1.6-millimeters  copper  graded  liner.  The  tin  and  copper  liner 
prevents  interference  by  lead  x-rays  [5]. 

The  detector  used  at  the  OSU  research  reactor  was  also  an  HPGe.  The  volume  of 
the  crystal  is  approximately  60  cubic  centimeters  [26]. 


Gamma  Spectroscopy 

Energy  and  efficiency  calibration  of  the  detector  to  a  known  source  is  essential  for 
accurate  measurements.  Once  calibrated,  an  unknown  source  can  be  measured  with  a 
high  degree  of  accuracy  provided  the  factors  that  can  impact  the  count  are  accounted  for. 

Absolute  efficiency  is  the  ratio  of  the  number  of  gamma  rays  counted  to  the 
number  of  gamma  rays  emitted  by  the  source.  It  is  represented  by 


Nc 

Nr 


fyfjjbs  ^int  ? 


(7) 


where  N c  is  the  number  of  gamma  rays  counted,  N E  is  the  number  of  gamma  rays 


emitted  by  the  source,  fy  is  the  self-absorption  correction  factor,  fw  is  the  detector 
window  correction  factor,  f  is  the  geometry  correction  factor,  fbs  is  the  back  scatter 


correction  factor  and  £int  is  the  intrinsic  efficiency  of  the  detector  itself. 
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Self-absorption  is  the  absorption  of  gamma  rays  by  the  source  itself.  The  mean 
free  path  of  a  gamma  in  a  material  is  the  average  distance  it  travels  before  interacting. 
Self-absorption  can  be  a  problem  if  the  mean  free  path  of  the  gamma  is  smaller  than  the 
thickness  of  the  foil.  The  mean  free  path  is  calculated  by 


X  = 


(8) 


where  A,  is  the  mean  free  path,  p  is  the  macroscopic  cross-section  and  p  is  the  density 
[28].  Self-absorption  is  usually  accounted  for  by  empirical  fonnulas  based  on 
experimental  results.  Equation  9  is  an  example  of  an  empirical  gamma  self-absorption 
correction  factor  for  gold  where  t  is  in  mils  [25]. 


4=i.o+ 


t 

239.5 


(9) 


The  detector  window  correction  factor ,fw,  corrects  for  the  number  of  gamma  rays 
absorbed  or  scattered  by  the  detector  window.  This  factor  is  normally  very  close  to  unity 
for  a  thin  window  and  the  detector  efficiency  typically  includes  this  correction  factor. 

Backscattering  of  gamma  rays  from  surrounding  materials  can  increase  the  count 
rate.  This  factor  ( fhs  )  is  also  normally  very  close  to  unity. 

The  counter  geometry  correction  factor,^,  corrects  for  the  fraction  of  the 
radiation  that  is  not  subtended  by  the  detector  and  is  defined  as 


4 


4;r 


(10) 
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where  Q  is  the  solid  angle  that  is  subtended  by  the  detector.  For  a  point  source  located 
along  the  axis  of  a  right  circular  cylindrical  detector,  Figure  2,  the  solid  angle  (in 
steradians)  is  given  by 

n  =  2*(i— =£=),  (ii) 

V d  +  a 

where  d  is  the  source-detector  distance  and  a  is  the  detector  radius. 


Figure  2.  Solid  angle  geometry  for  a  point  source. 


Figure  3.  Solid  angle  geometry  for  a  circular  disk. 

A  circular  disk  source  as  shown  in  Figure  3  has  a  more  complex  solution  to  the 
solid  angle.  Using  5  as  the  radius  of  the  source,  it  has  been  shown  that 
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(12) 


Q  = 


2  m'ee-dkJl(sk)Jl(ak) 

s  i  k 


where  the  Jx(x)  are  Bessel  functions  of  jc  [12]. 


Neutron  Sources 

Since  actual  nuclear  bomb  testing  is  unrealistic,  alternate  neutron  sources  must  be 
considered.  A  good  source  for  replicating  a  neutron  bomb  should  yield  as  many  fast 
neutrons  as  possible  and  be  capable  of  operating  in  a  controlled  outdoor  setting,  which 
would  be  ideal  to  represent  the  air-over-ground  transport  with  the  associated  ground  and 
sky  shine. 

Two  sources  were  used  in  this  thesis.  A  300  pCi  plutonium-beryllium  (Pu-Be) 
source  was  used  for  low  level  neutron  irradiation  and  the  500  kW  Ohio  State  University 
research  reactor  was  used  for  higher  level  neutron  irradiation.  Further  information  on  the 
Pu-Be  source  used  and  the  calculated  activity  is  contained  in  Appendix  A. 

The  OSU  research  reactor  is  an  open  pool  type  light  water  reactor.  The  fuel  is 
aluminum  clad  UjSf;  enriched  to  19.5%  U  with  a  '  U  loading  of  approximately 
3.9  kg.  The  pool  holds  approximately  5700  gallons  of  water  and  is  the  primary  coolant 
with  a  secondary  cooling  loop  containing  ethylene  glycol  and  water.  The  integral  neutron 
flux  near  the  core  has  been  measured  to  be  2  x  10  "  n  cm"  s"  at  full  power  [26]. 

Liquid/Air  Concept 

In  the  case  where  a  source  cannot  be  taken  outside,  the  distance  and  scattering 
mechanisms  must  be  simulated  or  engineered.  This  could  be  accomplished  by  placing  a 
depth  of  liquid  between  the  source  and  the  foil  packets.  The  theory  behind  this  is  that  the 
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liquid  is  much  denser  than  air  so  with  the  correct  geometry  the  number  of  molecules 
encountered  by  the  neutrons  passing  through  the  liquid  would  be  equivalent  to  the 
number  of  molecules  encountered  by  the  neutrons  passing  through  air  for  a  much  greater 
distance.  For  example,  the  density  of  dry  air  is  approximately  0.001 1  g  cm  3.  The 
density  of  liquid  air  is  approximately  0.89  g  cm"  .  The  ratio  of  the  density  of  liquid  air  to 
the  density  of  dry  air  yields  a  factor  of  809.  Thus,  1  centimeter  of  liquid  air  should 
appear  as  8.09  meters  of  dry  air  to  a  neutron.  Another  method  is  to  use  the  number  of 
atoms  in  a  cm  .  This  can  be  calculated  by 


PNAV(  o 
AW 


where  p  is  the  density,  NA  is  Avogadro’s  constant,  V  is  volume,  AW  is  the  atomic 
weight  and  (%)  is  the  percentage  of  occurrence  in  the  molecule.  Using  Equation  13, 
liquid  N20  has  6.3  x  10“  atoms  and  air,  modeled  as  78%  nitrogen  and  22%  oxygen  with 
a  density  of  1.1  mg/cm  ,  has  4.6  x  10  atoms  per  cm  .  This  gives  a  ratio  of  1371  or  1  cm 
of  liquid  N20  should  appear  as  13.71  meters  of  air.  Simulation  with  SCALE  indicated 
that  neither  method  is  correct  as  they  do  not  take  into  account  the  true  path  of  a  neutron 
moving  through  a  substance. 

Sky  shine,  the  scattering  of  neutrons  from  the  air  above,  should  still  be 
appropriately  replicated  in  liquid  air,  provided  the  height  of  the  liquid  air  is  sufficient  to 
reflect  all  neutrons  before  they  escape  from  the  container.  Other  liquids  closely 
resembling  air  theoretically  could  be  used  but  scattering  will  differ  depending  upon  the 
atomic  concentrations  (e.g.  heavier  atoms  scatter  differently  than  lighter  atoms).  Using 
liquid  to  replicate  air  in  depth  also  has  the  added  advantage  of  having  more  neutrons 
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available  at  the  experiment  if  the  same  source  is  used.  The  flux  from  an  isotropic  source 
decreases  by  the  square  of  the  distance  away  from  the  source.  Using  a  volume  of  liquid 
to  replicate  a  volume  of  air  removes  the  distance-squared  factor  because  of  the  smaller 
distance  traveled  by  the  neutron.  Thus,  a  much  higher  flux  can  be  attained  while  still 
accounting  for  scattering  through  a  large  distance  of  air.  To  apply  the  distance-squared 
factor  to  a  flux  measured  at  a  detector  near  the  source,  multiply  by 


d. 


d. 


2  ’ 


(14) 


where  dl  is  the  distance  from  the  source  to  the  near  detector  and  d2  is  the  distance  from 


the  source  to  the  point  desired. 


ORNL  Methodology 

ORNL  has  developed  a  computer  code  to  compare  an  unknown  neutron  spectrum 
to  a  database  of  known  neutron  weapon  spectra  and  rank-order  the  best  fit.  Prior  to  any 
detonation,  air-over- ground  adjoint  radiation  transport  calculations  are  performed  for 
each  type  of  known  foil  and  the  results  stored.  After  a  detonation,  as  many  foils  as 
possible  are  recovered  and  then  measured  to  determine  saturation  activity.  Once  the 
detonation  location  is  known,  the  calculated  saturated  activity  for  each  foil  is  obtained  by 
folding  all  the  known  weapon  leakage  spectra  with  the  pre-calculated  adjoint  fluences  at 
the  detonation  location.  The  calculated  and  measured  activities  can  now  be  compared 
and  the  most  probable  weapon  spectra  determined  by  least  squares  fitting  [22], 
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III.  Foil  Package  Design 


This  chapter  presents  the  design  of  the  foil  packages  starting  with  the  selection  of 
the  foil  materials  to  be  considered.  Simulations  using  SCALE  were  used  to  check  the 
effect  of  wire  geometry  on  activation  and  to  detennine  the  size  of  the  foils.  The  chapter 
concludes  with  the  packaging  of  the  selected  foils. 

Foil  Materials 

The  foils  considered  for  use  in  the  foil  packages  are  listed  in  Table  1.  It  is  based 
on  simulations  at  ORNL  and  includes  reactions  of  interest  that  give  as  complete  coverage 
of  the  neutron  spectrum  as  is  possible  [22]. 

Table  1.  Possible  candidate  materials  for  the  AFIDS  foil  package  [3,  15,  18,  28]. 


Reaction 

Flalf-Life 

Y  energy  TkeVl 

Threshold  IMeVl 

Wire  Form 

y  mean  free  nath  rcml 

27Al(n,a)24Na 

14.95  h 

1369,2754 

6.5 

Y 

7.0 

1 97  Au(n,2n) 1 96Au 

6.18  d 

356,  333,  426 

8.06 

Y 

0.16 

59Co(n,y)60Co 

5.271  y 

1333,  1173 

thermal 

Y 

2.0 

133Cs(n,p)133Xe 

5.243  d 

81 

3.0 

N 

0.14 

127I(n,2n)126I 

13.0  d 

389 

10.0 

N 

1.6 

115In(n,n')115In 

4.486  h 

497 

.369 

Y 

1.5 

139La(n,y)140La 

1.678  d 

1596,487,  816,  329 

thermal 

N 

0.92 

55Mn(n,2n)54Mn 

312.1  d 

835 

10.4 

Y 

2.1 

181Ta(n,y)182Ta 

114.43  d 

68,  1121 

thermal 

Y 

0.0052 

51V(n,a)48Sc 

43.7  h 

984.  1312,  1038 

9.0 

Y 

2.8 

58Ni(n,2n)57Ni 

35.6  h 

1378 

12.4 

Y 

2.1 

58Ni(n,p)58Co 

70.88  d 

811 

2.09 

Y 

1.6 

56Fe(n,p)56Mn 

2.578  h 

847.  1811,2113 

4.9 

Y 

1.9 

54Fe(n,p)54Mn 

312.1  d 

835 

thermal 

Y 

1.9 

63Cu(n,a)60Co 

5.271  y 

1333,  1173 

6.13 

Y 

3.5 

63Cu(n,y)64Cu 

12.701  h 

1346 

thermal 

Y 

3.5 
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The  mean  free  path  in  Table  1  is  calculated  using  Equation  8  and  the 
corresponding  underlined  gamma  energy.  Self-absorption  was  considered  unacceptable 
if  the  gamma  mean  free  path  was  shorter  than  the  foil  (wire)  diameter.  The  mean  free 
path  of  the  68  keV  gamma  from  tantalum  (branching  ratio  of  2.65)  is  too  short  but  the 
1.121  MeV  gamma  (branching  ratio  of  34.9)  has  a  mean  free  path  longer  than  1  mm. 

The  physical  characteristics  of  the  foil  are  a  significant  consideration  in  the  choice 
of  foils.  For  example,  cesium  with  a  3  MeV  threshold  reaction  is  not  available  as  a  wire 
and  is  pyrophoric.  Additionally,  cesium  activates  and  becomes  xenon,  a  noble  gas  that  is 
difficult  to  contain.  Iodine,  lanthanum  and  manganese  are  also  not  available  as  wires  but 
manganese  is  available  as  80%  manganese  /  20%  copper  wire  and  lanthanum  is  available 
in  a  foil  that  can  be  sliced  into  wires.  Iodine  was  not  chosen  for  this  design  but  could  be 
pelletized  and  considered  in  future  studies. 

Simulation  of  the  Effect  of  Wire  Shape  on  Neutron  Activation 

The  neutron  flux  is  attenuated  by  atoms  in  the  materials,  resulting  in  non-unifonn 
activation  of  the  sample.  For  this  reason,  foils  used  for  activation  analysis  are  typically 
very  thin  wires  or  circular  discs  that  are  oriented  toward  the  source  so  the  flux  penetrates 
through  the  material  with  little  attenuation  and  the  activation  is  relatively  unifonn.  Using 
the  point  detectors  available  in  SCAFE  within  a  wire  allows  comparison  of  the  activation. 
A  1  mm  thick  aluminum  wire  1  cm  long  was  simulated  500  m  from  a  14  MeV  point 
source  as  shown  in  Figure  4.  Nine  detectors  were  positioned  within  the  cross-section  of 
the  wire  as  shown  in  Figure  5. 
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Figure  4.  Foil  activation  simulation  geometry. 


SCALE  uses  the  Monte  Carlo  method  for  neutron  transport  and  thus  has  statistical 
error.  One  hundred  batches  of  1000  particles  were  required  to  keep  the  statistical  error 
below  10%.  If  the  variation  among  the  detectors  is  greater  than  the  statistical  error,  the 
flux  cannot  be  considered  uniform.  If  the  flux  is  uniform,  one  detector  can  be  placed 
anywhere  in  the  wire  for  future  simulations. 

The  total  response  column  in  Table  2  represents  the  activation  of  the  aluminum  at 
the  detector  location  as  a  result  of  the  flux  being  multiplied  by  the  cross-section  of 
aluminum.  The  fractional  standard  deviations  (fsd)  are  all  within  4.5%.  The  difference 
between  detector  7,  the  lowest,  and  detector  4,  the  highest,  is  just  1 .7%,  well  within  the 
fractional  standard  deviation.  This  allows  one  detector  to  be  placed  in  the  wire  for  future 
simulations. 


21 


Table  2.  Partial  SCALE  output  of  1mm  A1  wire  at  500  meters. 

responses(detector) 


detector 

uncoil 

fsd 

total 

fsd 

response 

uncoil 

response 

total 

1 

9.6  x  1010 

1.8  x  10-3 

6.5  x  1011 

4.4  x  10'2 

2 

9.2  x  1010 

1.9  x  10-3 

6.5  x  1011 

4.4  x  10'2 

3 

9.3  x  1010 

2.0  x  10-3 

6.5  x  1011 

4.5  x  10'2 

4 

9.8  x  1010 

1.8  x  10-3 

6.6  x  1011 

4.5  x  10'2 

5 

9.4  x  1010 

1.8  x  10-3 

6.5  x  1011 

4.4  x  10'2 

6 

9.3  x  1010 

2.0  x  10-3 

6.5  x  1011 

4.4  x  10'2 

7 

9.3  x  1010 

2.0  x  10-3 

6.5  x  1011 

4.3  x  10'2 

8 

9.3  x  1010 

2.0  x  10'3 

6.5  x  1011 

4.3  x  10'2 

9 

9.3  x  1010 

2.0  x  10'3 

6.5  x  1011 

4.4  x  10'2 

Foil  Size 

Another  limitation  on  the  foils  was  that  activity  had  to  be  measurable  four  days 
after  activation.  Simulations  were  used  to  determine  the  size  of  foil  needed  to  obtain  this. 
The  simulations  were  run  on  each  type  of  foil  at  500  m  from  the  source.  Air  was 
configured  as  78%  nitrogen  and  22%  oxygen  with  a  density  of  1 . 1  mg/cm  .  The  foils 
were  right  circular  cylinders,  1  mm  in  diameter,  1  cm  long  and  centered  at  the  simulated 
distance  standing  upright  as  shown  in  Figure  4.  A  point  source  was  simulated  having 
10“  neutrons  and  energy  of  14  MeV. 

The  SCALE  results  were  used  to  calculate  the  activity  four  days  after  the 
activation.  The  sizes  of  the  wires  were  adjusted  to  attain  an  activity  between  the 
minimum  for  the  particular  wire  and  1000  Bq.  The  complete  spreadsheet  is  in  Appendix 
D  and  the  significant  results  are  presented  in  Table  3. 
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Table  3.  Determining  foil  size  from  projected  activity. 


(ct<)>)  SCALE 

radius 

length 

mass 

Initial  A 

VA 

A  at  +4  days 

Min.  A 

Al 

500  m 

[n/s] 

[cm] 

[cm] 

[mg] 

[Bq] 

[days] 

[Bq] 

[Bq] 

6.26  x  1 0'13 

0.05 

1.00 

21.20 

3815.97 

0.62 

44.56 

10.0 

330  m 

3.24  x  10"12 

0.05 

0.05 

1.00 

932.27 

0.62 

10.89 

10.0 

700  m 

2.07  x  1 0'13 

0.05 

1.00 

21.20 

1259.02 

0.62 

14.70 

10.0 

800  m 

7.64  x  1 0'14 

0.05 

1.00 

21.20 

465.41 

0.62 

5.44 

10.0 

Au 

500  m 

9.46  x  10"12 

0.05 

0.30 

45.47 

1707.56 

6.18 

1090.37 

12.4 

1100  m 

1.85  x  1 0'13 

0.05 

0.30 

45.47 

33.33 

6.18 

21.28 

12.4 

1200  m 

6.22  x  1 0'14 

0.05 

0.30 

45.47 

11.22 

6.18 

7.17 

12.4 

Co 

500  m 

6.24  x  1 0'11 

0.05 

1.00 

69.90 

185.90 

1923.92 

185.63 

10.0 

1100  m 

3.80  x  1 0'12 

0.05 

1.00 

69.90 

11.33 

1923.92 

11.31 

10.0 

1200  m 

In 

100  m 

8.38  x  1 0'13 

0.05 

1.00 

69.90 

2.50 

1923.92 

2.49 

10.0 

1.22  x  1 0'9 

0.05 

1.00 

57.41 

1.58  x  107 

0.19 

5.73 

23.0 

500  m 

2.49  x  1 0'11 

0.05 

1.00 

57.41 

3.22  x  106 

0.19 

0.12 

23.0 

La 

500  m 

1.15  x  10"11 

0.05 

0.50 

24.13 

5747.28 

1.68 

1101.60 

10.5 

1600  m 

2.72  x  1 0'13 

0.05 

0.50 

24.13 

136.26 

1.68 

26.12 

10.5 

1700  m 

1.32  x  1 0'14 

0.05 

0.50 

24.13 

6.62 

1.68 

1.27 

10.5 

Mn  80% 

500  m 

2.71  x  10'12 

0.05 

1.50 

68.80 

52.60 

312.10 

52.14 

10.0 

700  m 

5.23  x  1 0'13 

0.05 

1.50 

68.80 

10.14 

312.10 

10.05 

10.0 

Ta 

500  m 

3.83  x  1 0'10 

0.02 

0.50 

9.44 

844.82 

114.43 

824.60 

28.7 

1000  m 

6.65  x  10'11 

0.02 

0.50 

9.44 

146.42 

114.43 

142.91 

28.7 

1100  m 

V 

500  m 

1.33  x  1 0'11 

0.02 

0.50 

9.44 

29.33 

114.43 

28.63 

28.7 

6.23  x  1014 

0.05 

1.00 

47.99 

155.61 

1.82 

33.95 

10.0 

600  m 

3.38  x  1014 

0.05 

1.00 

47.99 

84.41 

1.82 

18.42 

10.0 

700  m 

1.60  x  1 0'14 

0.05 

1.00 

47.99 

40.00 

1.82 

8.73 

10.0 

Ni 

500  m 

cross-section  for  (n,p) 
3.67x1  O'12  0.05  1.00 

69.92 

298.12 

70.88 

286.68 

10.1 

1000  m 

1.57  x  1 0'13 

0.05 

1.00 

69.92 

12.74 

70.88 

12.25 

10.1 

1100  m 

9.82  x  1 0'14 

0.05 

1.00 

69.92 

7.97 

70.88 

7.67 

10.1 

Cu 

250  m 

cross-section  for  (n,a) 

1.61  xIO-12  0.05  2.00 

140.74 

8.93 

1923.92 

8.92 

10.0 

500  m 

1.51  x  IQ'13 

0.05 

2.00 

140.74 

0.84 

1923.92 

0.84 

10.0 

The  wires  chosen  for  further  experimentation  are  listed  in  Table  4  along  with  their 
dimensions.  Iron  did  not  have  enough  activation  with  reasonable  size  to  be  measured 
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four  days  after  an  event  so  was  not  pursued  as  a  foil  in  further  experiments.  The  activity 
of  indium  was  also  below  the  minimum  of  23  Bq  due  to  the  short  half-life  but  could  still 
be  useful  if  collected  early  enough.  Most  of  the  wires  were  chosen  with  a  length  of  1  cm 
but  gold  was  shortened  to  3  mm  to  reduce  its  activation  to  near  1000  Bq,  this  also  kept 
the  cost  down.  The  foil  length  was  limited  to  3  mm  to  keep  it  upright  in  the  glass  tubes 
used  to  package  the  foils.  The  length  the  manganese/copper  wire  was  set  at  1.5  cm  to 
assist  the  activity  of  both.  The  tantalum  wire  was  set  to  0.5  cm  to  keep  activity  under 
1000  Bq.  Copper  did  not  meet  the  minimum  activation  but  two  wires  of  1  cm  each  were 
chosen  as  it  is  relatively  inexpensive.  Table  3  also  contains  the  expected  distance  from 
the  source  the  foil  was  expected  to  perform.  For  example,  vanadium  should  perform  as 
designed  to  between  600  and  700  meters  as  it  activated  to  18.42  Bq  at  600  meters  but 
only  8.73  Bq  at  700  meters. 


Table  4.  Foils  selected  for  experiment. 


diameter 

length  each 

Element 

rmml 

Tcml 

A1 

1 

1 

Au 

1 

0.3 

Co 

1 

1 

In 

1 

1 

La 

1 

0.5 

Mn  -  Cu 

1 

1.5 

Ta 

0.38 

0.5 

V 

1 

1 

Ni 

1 

1 

Cu 

1 

2 

A  check  on  the  simulation  was  done  using  MCNP  to  simulate  a  1  mg  A1  wire  at 
500  meters.  The  MCNP  results  and  corresponding  calculations  are  in  Appendix  D. 
SCALE  gave  1.4  x  107  atoms  activated  and  MCNP  yielded  1.7  x  107.  The  results  were 
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within  18%.  This  check  was  done  to  see  if  SCALE  was  being  implemented  properly  and 
results  were  realistic  before  ordering  the  wires. 

Foil  Packaging 

The  geometry  shown  in  Figure  6  meets  the  360-degree  coverage  limitation  while 
keeping  the  foils  from  shielding  themselves  or  each  other.  Initial  planning  has  the  foil 
packages  placed  on  high  structures  such  as  building  tops  and  cell  phone  towers  to 
maximize  exposure. 


Figure  6.  Foil  package. 

The  packaging  should  also  not  shield  the  foils  or  it  may  change  the  spectrum  for 
activation.  Silicon  carbide  would  make  an  ideal  packaging  material  with  a  melting  point 
of  1825°C,  a  modulus  of  rupture  of  9500  psi,  and  a  low  absorption  cross-section  [24], 

The  glass  tubes  used  in  this  thesis  contained  sodium  that  activated  and  were  too  fragile  to 
survive  in  the  field. 
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IV.  Foil  Package  Activation  and  Spectroscopy 


This  chapter  presents  the  results  of  foil  package  activation  and  spectroscopy 
measurements.  The  gamma  spectroscopy  systems  used  were  calibrated  against  each 
other  and  the  accuracy  of  measuring  wires  as  point  sources  investigated.  The  foils  were 
activated  to  a  low  level  to  check  for  complications  in  gamma  spectroscopy  and  then  the 
foil  packages  were  activated  to  a  higher  level  to  test  performance. 

Gamma  Spectroscopy 

Foils  had  to  be  activated  to  perform  gamma  spectroscopy  experiments.  They 
were  activated  locally  using  a  Pu-Be  source.  The  Pu-Be  source  is  contained  in  a  barrel 
approximately  45  cm  high  and  36  cm  in  diameter.  The  barrel  is  filled  with  paraffin  and 
has  a  hollow  center  core  approximately  10  cm  in  diameter.  The  source  is  stored  at  the 
bottom  of  the  hole,  approximately  21.5  cm  from  the  end  of  the  barrel.  Wires  were 
activated  with  the  Pu-Be  source  container  lying  on  its  side  and  the  wires  contained  in  foil 
packets  located  over  the  open  end  as  shown  in  Figure  7.  The  neutron  flux  at  the  foil 
location  was  1.9  x  10'  n  cm''  s'  as  calculated  in  Appendix  A. 
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Spectroscopy  System  Calibration. 

Accurate  spectroscopy  requires  energy  and  efficiency  calibration.  If  the 
calibration  is  done  in  the  same  geometry  as  the  measuring  and  the  sources  are  similar,  no 
correction  factors  need  to  be  applied.  The  HPGe  detector  system  with  the  Genie  2000 
software  [6]  was  calibrated  using  a  multi-nuclide  source  consisting  of  the  isotopes  in 
Table  5  with  a  reference  date  of  15  July  2004.  The  gamma  spectrum  was  measured  using 
the  same  geometry  proposed  for  the  wires  (1  mm  off  the  detector  face),  including 
background  subtraction. 

Table  5.  Multinuclide  source  for  calibration. 


Nuclide 

Half  Life  Id) 

Initial  Activity  (v/si 

Y  (keV) 

24 'Am 

157742.1 

387.9 

60.0 

109Cd 

462.2 

572.2 

88.0 

57Co 

271.8 

484.3 

122.1 

I39Ce 

137.6 

580.3 

165.9 

203Hg 

46.6 

1845 

279.2 

113Sn 

115.1 

1797 

391.7 

85  Sr 

64.8 

3431 

514.0 

I37Cs 

11012.1 

2165 

661.6 

>- 

oo 

oo 

106.6 

5151 

898.0 

60Co 

1924.3 

2989 

1173.2 

60Co 

1924.3 

2993 

1332.5 

88y 

106.6 

5447 

1836.0 

A  1  cm  Cu-Mn  wire  (1  mm  diameter)  was  activated  for  30  hours  to  compare 
AFIT’s  spectroscopy  system  with  OSU’s  spectroscopy  system.  The  wire  was  chosen  for 
the  large  thennal  cross-section  of  the  55Mn(n,y)56Mn  reaction  and  the  2.578  hour  half  life. 
This  combination  allowed  the  wire  to  be  irradiated  to  saturation  with  the  Pu-Be  source 
locally  and  then  transported  to  OSU  for  gamma  spectroscopy  with  measurable  activity. 
The  OSU  system  had  been  calibrated  with  a  1  mm  diameter,  1  cm  long  cobalt  wire  source 
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from  the  National  Institute  of  Standards  and  Technology  (NIST)  [26],  The  activity 
calculated  by  the  OSU  spectroscopy  system  for  the  Cu-Mn  wire  was  5.53  becquerel  (Bq) 
±  15.95%  at  the  time  the  wire  was  removed  from  the  Pu-Be  source.  After  another 
activation  to  saturation  with  the  Pu-Be  source,  the  wire  was  measured  on  the  AFIT 
spectroscopy  system.  The  activity  of  4.775  ±  0.3 13  Bq  was  within  the  15.95%  error  of 
the  5.53  Bq  measured  at  OSU.  The  experiment  was  repeated  using  a  gold  wire  resulting 
in  activities  of  4916  Bq  ±  3.25%  at  OSU  and  5 1 16  ±  143  Bq  at  AFIT.  These  results  were 
also  within  the  error  limits.  The  efficiency  calibration  used  is  shown  in  Figure  8. 


( Efficiency  Calibration  Curves 


Energy  (keV) 


Datasource:  C:\GENIE2K\CAMFILES\KimmelVAu-0SU.CNF 
In(Eff)  =  -2.442e+000  +  6.1 1 2e-001  Xx  - 1 ,288e-001  Xxrt2  •  4.220e-001  xxA3 
+  4.378e-001  xxA4  - 1 .1 59e-001  Xx"5 
where:  x  =  ln(  9.790e+002/E) 


OK 
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Figure  8.  GENIE  2000  output  showing  the  calibration  equation  fit 
from  the  multinuclide  source. 


Wire  vs.  Point. 

A  comparison  of  the  solid  angle  using  Equations  1 1  and  12  (Appendix  C) 
indicated  that  a  wire  source  could  not  be  measured  exactly  as  a  point  source.  However,  at 
1  mm  off  the  detector  surface,  the  difference  in  solid  angle  between  a  point  and  a  1  cm 
wire  is  less  than  0.1%.  Since  this  amount  of  error  is  negligible  when  compared  to  the 
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normal  statistical  error  of  the  spectroscopy  system  (greater  than  1%),  the  wire  could  be 
calibrated  and  measured  as  though  it  were  a  point  source. 

To  experimentally  validate  this  conclusion,  a  1  cm  gold  wire  was  activated  by  the 
Pu-Be  source  and  measured.  The  activity  was  measured  at  5.184  x  10'4  ±  1.84  x  10'5  pCi. 
The  wire  was  then  rolled  into  a  sphere  approximately  2  mm  in  diameter  to  more  closely 
resemble  a  point  source  and  measured  again.  The  result  was  an  activity  of  4.563  x  10’4  ± 
1.58  x  10"5  pCi.  These  results  do  not  correspond  within  error  bounds  but  when 
Equation  9  is  applied  to  compensate  for  self-absorption  in  gold,  the  results  become 
5.321  x  10  4  ±  1.84  x  10°  pCi.  These  results  correspond  within  the  error  and  confirm  that 
the  difference  between  measuring  a  point  source  and  a  1  cm  wire  source  at  the  face  of  the 
detector  is  insignificant,  backing  up  the  solid  angle  comparison  in  Appendix  C. 

Low  Level  Activation  Measurements 

The  foils  were  activated  with  a  15-hour  irradiation  by  the  Pu-Be  source  to  check 
for  gamma  producing  reactions  that  could  interfere  with  the  gamma  spectroscopy  of  the 
reactions  of  interest.  The  foils  received  a  neutron  fluence  on  the  order  of  1  x  10  n  cm" 
and  only  three  reactions  were  detected  in  two  hour  counts  as  shown  in  Table  6.  The 
indium  reaction  (54  minute  half-life)  and  manganese  reaction  (2.5  hour  half-life)  would 
both  decay  prior  to  gamma  spectroscopy  on  the  reactions  of  interest.  The  gold  reaction 
(2.7  day  half-life)  would  still  be  present  at  foil  recovery  but  the  reaction  of  interest  for 
gold  (n,2n)  has  a  half-life  of  6.2  days  so  the  longer  the  delay  in  counting,  the  easier  it 
would  be  to  detect  the  (n,2n)  reaction  over  the  (n,y)  reaction. 
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Table  6.  Pu-Be  activation  results. 


Element 

Primary  Gammas  Detected  (keV) 

Reaction  Activity  ( uCi) 

Ni 

None 

La 

None 

Cu 

None 

V 

None 

Ta 

None 

Co 

None 

In 

417,  819,  1098,  1294,2113 

n,y  6.67  x  10'5 

A1 

None 

Au 

411 

n,y  4.53  x  10'6 

Mn 

846,  1810 

n,y  4.48  x  10  5 

The  four  most  common  neutron  reactions,  (n,y),  (n,2n),  (n,p)  and  (n,a),  were 
researched  for  each  foil  and  the  half-life,  decay  method  and  percentage,  branching  ratio 
and  gamma  energy  for  each  were  tabulated  in  Appendix  B.  The  list  was  then  sorted  by 
gamma  energy  and  possible  spectroscopy  conflicts  could  be  isolated.  The  screening 
criterion  was  that  any  two  gammas  with  peaks  within  2  keV  of  each  other  would  not  be 
measured  together  since  the  resolution  of  an  HPGe  detector  is  ~2  keV  FWHM  at 
1 .3  MeV.  Copper,  cobalt  and  nickel  have  reactions  that  produce  the  same  nuclide  and 
tantalum  and  gold  have  reactions  that  produce  several  gammas  that  are  within  2  keV  of 
each  other.  Nine  more  combinations  of  materials  had  gammas  within  2  keV  but  by 
grouping  the  wires  into  three  groups  for  measuring,  all  but  one  combination  can  be 
avoided.  The  136  keV  gamma  is  common  with  nickel  and  gold.  The  nickel  is  a  second- 
generation  decay  and  gamma  from  the  gold  comes  from  a  metastable  state  and  has  a  1.3% 
branching  ratio.  Both  reactions  produce  other  gammas  with  much  higher  branching 
ratios.  Measuring  each  group  of  wires  at  once  reduces  the  measurements  required  for  a 
foil  package  from  ten  to  three.  Group  1  of  the  wires  for  measuring  is  cobalt,  tantalum, 
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manganese  and  aluminum;  group  2  is  indium,  lanthanum  and  copper;  and  group  3  is 
nickel,  vanadium  and  gold. 

High  Level  Activation  Measurements 

To  verify  the  foils’  reactions  of  interest  would  activate  at  the  designed  levels  they 
were  irradiated  in  the  rabbit  tube  of  the  OSU  research  reactor.  The  rabbit  tube 
pneumatically  locates  the  item  for  irradiation  next  to  the  core.  The  operator  can  manually 
insert  and  eject  the  item  for  irradiation  or  a  timer  can  be  used  to  eject  automatically  at  the 
preset  time.  The  rabbit  tube’s  neutron  energy  spectrum  from  previous  experimentation  is 
shown  in  Figure  9. 


The  gold,  aluminum,  copper  and  manganese  foils  in  a  package  were  irradiated  in 
the  rabbit  tube  for  1000  seconds  at  20%  power  to  expose  them  to  a  fluence  of 
1.8  x  10  n  cm'  for  fast  and  epithermal  neutrons  (over  0.5  eV).  The  indium,  vanadium 

and  nickel  foils  in  a  package  were  irradiated  in  the  rabbit  tube  for  1000  seconds  at  10% 
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13  2 

power  to  expose  them  to  a  fluence  of  9  x  10  n  cm'".  The  foil  packages  were  wrapped  in 
cadmium  to  partially  inhibit  thennal  activation  so  the  packages  could  be  safely  handled 
when  removing  them  from  the  rabbit  tube.  The  results  are  tabulated  in  Table  7.  The 
(n,a)  reaction  of  copper  the  (n,2n)  reaction  of  manganese  were  not  detected.  These  both 
have  relatively  long  half-lives  of  5.27  years  and  312  days  respectively.  Gold  was  initially 
measured  two  days  after  irradiation  and  had  so  much  activity  that  the  dead  time  was  at 
98%.  This  was  a  product  of  the  (n,y)  reaction  and  not  the  reaction  of  interest  (n,2n).  The 
dead  time  was  reduced  to  25%  by  moving  the  foil  13  cm  from  the  detector.  Another 
count  was  made  19  days  after  irradiation  and  the  dead  time  was  down  to  7%  with  the 
wire  placed  1  mm  from  the  face  of  the  detector.  The  peak  of  interest  (n,2n)  was  also 
relatively  larger  than  the  shorter  lived  (n,y)  reaction. 

Table  7.  Rabbit  tube  activation  results. 


Element 

ROl 

T!4  (h) 

Delay  Time  (h) 

Measured  Act  (pCi) 

Act  at 

Activation  (|iCi) 

Au 

n,2n 

148 

475.77 

3.04  x  10'5 

2.82  x 

io-4 

Au 

n,Y 

64.685 

475.77 

1.35  x  10'1 

2.21  x 

101 

A1 

n,a 

14.95 

43.05 

2.45  x  10'7 

1.80  x 

10'6 

Cu 

n,a 

46205 

43.05 

N/A 

N/A 

Cu 

n,y 

12.701 

43.05 

1.12 

1.18  x 

101 

Mn 

n,2n 

7490 

43.05 

N/A 

N/A 

Mn 

n,y 

2.578 

43.05 

5.43  x  10'4 

5.78  x 

101 

In 

n,n' 

4.485 

18.133 

1.76  x  10'1 

2.90 

In 

n,Y 

0.903 

18.133 

4.40  x  10'4 

4.88  x 

102 

V 

n,a 

43.7 

18.133 

3.17  x  10"5 

4.23  x 

10'5 

Ni 

n,2n 

35.6 

18.133 

9.55  x  10'6 

1.36  x 

10'5 

Ni 

n,p 

1701.1 

18.133 

7.60  x  10'3 

7.66  x 

10'3 
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V.  Liquid  Air  Experiments 


This  chapter  investigates  the  concept  of  using  liquid  to  replicate  large  volumes  of 
air.  The  investigation  consisted  of  a  simulation  of  the  general  concept,  the  design  of  a 
physical  experiment,  a  simulation  of  the  experiment,  and  the  conduct  of  the  experiment. 


Liquid/Air  Concept  Simulation 

The  concept  of  using  compressed  gasses  or  liquids  to  represent  distances  of  air 
was  investigated  using  SCALE  to  simulate  neutron  transport  through  7  cm  of  liquid  N2O 
as  shown  in  Figure  10.  A  point  source  was  modeled  as  monoenergetic  (14  MeV)  with  a 
source  magnification  of  10.  All  boundaries  were  modeled  as  spheres  so  nothing  would 
affect  the  transport  of  neutrons  through  air  except  the  liquid  N2O.  A  point  detector  was 
placed  240  cm  from  the  source  to  measure  the  neutron  flux  after  transporting  through  the 
liquid  N20.  To  get  neutron  flux  as  the  output  of  SCALE,  the  response  function  was  set  to 
one.  The  N20  was  modeled  as  67%  nitrogen,  33%  oxygen  with  density  of 

•3 

1 .53  g  cm'  [31].  The  air  was  modeled  as  78%  nitrogen  and  22%  oxygen  with  a  density 
of  1 . 1  mg  cm'  . 


Not  to  Scale 


Figure  10.  Nitrous  oxide  simulation. 


33 


The  simulation  was  repeated  by  replacing  the  N20  with  liquid  air  modeled  the 
same  as  the  dry  air  but  with  density  of  1 .53  g  cm'  .  The  simulation  was  again  repeated 
replacing  the  N20  with  the  modeled  air  and  placing  point  detectors  between  93  and  175 
meters  to  measure  the  neutron  flux  after  transporting  through  air. 

Using  the  number  of  atoms  per  unit  volume,  7  cm  of  N20  should  be  equivalent  to 
96  meters  of  air.  The  results  of  the  simulations,  graphed  in  Figure  1 1  with  the  dispersion 
factor  (Equation  14)  applied  to  the  N20  and  liquid  air  results,  lead  us  to  believe  that  it  is 
not  quite  that  simple.  If  the  7  cm  of  N20  was  equivalent  to  96  meters  of  air  for  neutron 
transport,  the  results  should  have  crossed  at  96  meters  on  the  graph. 


Nitrous  vs.  Air  vs.  Liquid  Air 


Figure  1 1 .  Comparison  of  SCALE  simulations  for  N20,  air  and  liquid  air. 

Liquid/Air  Experiment  Design 

An  experiment  using  the  beam  port  of  the  OSU  research  reactor  was  designed 
with  the  original  goal  of  sending  foil  activity  measurements  to  ORNL  by  using  N20  to 
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simulate  air  over  ground  neutron  transport  while  representing  different  volumes  of  air 
(e.g.  19  cm  of  N2O  in  the  cylinder  represents  264  meters  of  air).  Experiments  and 
simulations  were  carried  out  concurrently  and  the  results  used  to  shape  further 
experiments  and  simulations.  The  results  of  the  liquid/air  concept  simulation  changed  the 
original  goal  and  experiment  to  comparing  activation  results  using  N?0  to  results  without 
N?0  to  detennine  an  empirical  fit.  The  foil  activations  would  be  accomplished  in  two 
consecutive  reactor  runs.  The  first  run  would  be  without  N2O  with  the  foils  located 
21.5  cm  inline  and  outside  the  end  of  the  beam  port.  For  the  second  run,  the  N2O 
cylinder  would  be  attached  to  the  end  of  the  beam  port  and  the  foils  located  1.5  cm 
beyond  the  cylinder  and  inline  with  the  beam  port. 

OSU’s  research  reactor  has  the  capability  of  producing  a  wide  spectrum  of 
neutrons  as  shown  in  Figure  9.  The  integral  neutron  flux  near  the  core  has  been 
measured  to  be  2  x  10  "  n  cm'"  s'  and  10  n  cm'  s'  at  the  outside  end  of  the  beam  port  at 
full  power  [26].  The  reactor  also  has  strict  operating  restrictions  on  the  amount  of 
radiation  (2  mrem  h'1)  that  can  escape  the  facility.  The  open  beam  port  is  an  avenue  for 
radiation  that  must  be  shielded  to  maintain  safe  levels  of  exposure. 

Additionally,  to  replicate  a  nuclear  detonation  the  source  must  appear  as  a  point. 
Therefore,  two  polyethylene  plugs,  30  cm  each,  with  a  38  mm  hole  in  the  center  were 
added  to  the  beam  port  to  collimate  the  source  and  make  it  appear  as  a  point  as  shown  in 
Figure  12. 
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|  Foil  Packet 


Not  To  Scale 


Figure  12.  Collimated  OSU  source  experiment. 

The  plugs  were  inserted  82.4  cm  from  the  outside  of  the  beam  port  to  make  the 
source  appear  as  a  point  to  the  six  inch  opening  at  the  end  of  the  beam  port.  The  method 
of  using  similar  triangles  was  used  to  detennine  the  distance  of  insertion  required  as 
shown  in  Figure  13. 


142.42  cm 

◄ - ► 


15.24  cm 


Similar  Triangles 
1.9  7.12 

30  ~  X  +  30  M - - ► 

Not  to  Scale  30  cm  X 

Figure  13.  Collimating  the  OSU  reactor  to  appear  as  a  point  source. 

Cadmium  was  added  to  the  outside  end  of  the  plugs  and  the  beam  port  between 
the  end  of  the  plugs  and  the  outside  opening  was  lined  with  cadmium  to  prevent  thermal 
neutrons  from  scattering  back  in.  The  hydrogen  in  the  polyethylene  (CfU)  and  the 
Barytes  cement  containment  structure  slows  the  neutrons  down  by  elastic  scattering  and 
the  cadmium,  with  a  huge  thermal  absorption  cross-section,  captures  them. 
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Liquid/Air  Experiment  Simulation 

The  OSU  research  reactor  beam  port  experiment  was  simulated  using  SCALE. 
The  simulation  was  done  using  the  geometry  depicted  in  Figure  14.  Experimental 
components  were  then  added  as  needed  to  simulate  individual  experiments. 


Figure  14.  OSU  reactor  beam  port  simulation  geometry. 

The  source  was  modeled  as  a  48  cm  cube  surrounded  by  water.  The  inner  wall  of 
the  cement  containment  model  was  modeled  as  a  sphere  56.5  cm  from  the  center  of  the 
core.  The  outer  wall  was  modeled  as  a  sphere  230  cm  from  the  center  of  the  core.  The 
cement  was  modeled  with  a  density  of  3.3  g  cm'  corresponding  with  Barytes 
concrete  [28],  The  beam  port  was  modeled  as  a  right  circular  cylinder  consisting  of  two 
parts  with  different  diameters.  The  part  closest  to  the  core  started  1  cm  from  the  core, 
extended  159  cm  horizontally  and  had  a  radius  of  7.62  cm.  The  outside  portion  began 
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where  the  first  part  ended  (184  cm  from  center  of  core)  and  extended  to  the  outside  of  the 
cement  containment  structure  (46  cm).  The  air  in  the  beam  port  and  outside  the  cement 
containment  structure  was  modeled  as  78%  nitrogen  and  22%  oxygen  with  a  density  of 
1.1  mg/cm3.  Point  detectors  were  centered  in  the  beam  port  1  cm  from  the  end  closest  to 
the  core,  1  cm  outside  the  end  of  the  beam  port  and  21.5  cm  off  the  end  of  the  beam  port 
where  the  foils  would  be  placed.  The  response  function  was  set  to  one  so  that  the  output 
of  SCALE  would  be  the  flux.  The  source  spectrum  used  was  from  previous  experimental 
data  taken  in  the  central  irradiation  facility  (CIF)  located  in  the  center  of  the  core  [10]. 
The  source  multiplication  factor  was  3.36x  1016  and  was  determined  by  trial  and  error  to 
get  the  integral  flux  at  the  end  of  the  beam  port  closest  to  the  reactor  core  to 
3.6  x  10  “  n  cm'  s'  ±  19.5%,  corresponding  to  previous  experimental  data  taken  from  the 
rabbit  tube  [29]. 

The  neutron  flux  from  the  simulation  at  21.5  cm  outside  the  beam  port,  the  foil 
placement  point,  was  1.903  x  10  n  cm'-  s'  ±  17.2%.  Using  this  result  and  operating  at 
450  kW,  it  would  require  about  an  hour  and  a  half  to  reach  a  fluence  of  10  n  cm'  and 
over  14  hours  for  1014  n  cm"2.  This  simulation’s  resulting  neutron  flux  was  an  order  of 
magnitude  higher  than  what  had  been  previously  measured. 

Unsafe  levels  of  gamma  and  neutron  radiation  were  measured  in  the  reactor 
building  upon  experimentation  with  the  collimating  polyethylene  plugs.  Because  of  this, 
all  further  beam  port  simulations  replaced  the  collimating  polyethylene  plugs  with  a  solid 
polyethylene  plug  45  cm  long  in  the  outside  end  of  the  beam  port  to  replicate  actual 
experiments. 
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Keeping  all  parameters  the  same  and  exchanging  the  collimating  polyethylene 
plugs  with  the  solid  plug  at  the  end  of  the  beam  port  gave  a  simulated  neutron  flux  at  the 
foil  locations  of  1.4625  x  10  n  cm'  s'  ±  17.2%.  This  flux  was  five  orders  of  magnitude 
less  than  the  flux  from  the  simulation  with  the  collimating  polyethylene  plugs  and  should 
result  in  a  safe  level  of  radiation  in  the  reactor  building. 

The  effects  of  nitrous  oxide  were  simulated  using  an  N2O  container  simulated  by 
a  right  circular  cylinder  of  4130  steel  with  sidewalls  of  0.5  cm  [32].  The  cylinder  was 
20  cm  in  diameter  and  5 1  cm  high.  It  was  placed  perpendicular  against  the  end  of  the 
beam  port  so  the  outside  edge  was  250  cm  from  the  center  of  the  core.  The  detector  was 
placed  1 .5  cm  from  the  cylinder  centered  on  the  beam  port  or  25 1 .5  cm  from  the  center  of 
the  core  as  shown  in  Figure  15.  The  4130  steel  was  modeled  as  in  Table  8.  The  N2O  was 
modeled  as  67%  nitrogen,  33%  oxygen  with  density  of  1.53  g  cm" .  The  resulting  flux  at 
the  foil  location  was  2.2041  x  103  n  cm"2  s'1  ±  17.5%. 


|  Foil  Packet 


Figure  15.  OSU  beam  port  characterization  with  poly  plug  and  N20. 
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Table  8.  4130  steel  composition  [9]. 


Element 

Percentage 

Iron 

97.56 

Carbon 

0.30 

Manganese 

0.65 

Phosphorus 

0.04 

Sulfur 

0.05 

Silicon 

0.25 

Chromium 

0.95 

Molybdenum 

0.20 

This  flux  was  too  low  to  activate  any  of  the  reactions  of  interest  for  AFIDS  to 
measurable  amounts  but  should  be  enough  to  activate  the  (n,y)  reactions  of  gold  and 
copper  to  characterize  the  neutron  spectrum  if  the  escaping  radiation  could  be  controlled 
enough  to  allow  a  high  power  lengthy  run.  This  spectrum  could  then  be  compared  to  the 
spectrum  obtained  from  the  run  without  the  N2O  to  get  more  insight  on  the  effect  of  N2O 
on  neutron  transport. 

Experiment 

Characterization  of  neutron  spectra  is  done  at  the  OSU  research  reactor  using  the 
SAND-II  program.  The  SAND-II  program  determines  neutron  spectra  by  multiple  foil 
integration  using  the  iterative  method  and  generates  a  620-group  neutron  energy 
spectrum.  The  required  inputs  are  the  saturation  activity  per  nucleus  of  the 
characterization  foils  used.  The  characterization  foils  used  at  the  OSU  research  reactor 
consisted  of  gold,  cobalt  and  copper  in  two  sets.  One  set  was  wrapped  in  cadmium  and 
the  other  was  bare.  The  spreadsheet  in  Table  9  was  created  for  calculating  the  saturation 
activity  per  nucleus. 


40 


Table  9.  SAND-II  input  calculations. 


wire/foil 

*  mass 
(mg) 

%  AW 

abund.  (g/mol) 

Target 

Nuclei 

Isotope 

of 

Interest 

* 

Activity 

(dps) 

T!4  (m) 

X  (m'1) 

*  Time 
Irradiated 
(m) 

Saturated 

Activity 

(dps/nucleus) 

Bare 

Au 

8.99 

100 

196.97 

2.75  x  1019 

I98Au 

110.59 

3881.09 

0.00018 

15 

1.50  x  10'15 

Co 

1 

100 

58.93 

1.02  x  1019 

60Co 

100 

2772276  2.5  x  10-7 

15 

2.61  x  10'12 

Cu 

1 

69.17 

62.93 

9.57  x  1018 

64Cu 

100 

762.06 

0.00091 

15 

7.71  x  10'16 

Fe 

1 

0.282 

57.93 

1.04  x  1019 

59Fe 

100 

64080 

1.1  x  10-5 

15 

5.93  x  10'14 

Cd 

covered 

An 

1 

100 

196.97 

3.06  x  1018 

l98Au 

100 

3881.09 

0.00018 

15 

1.22  x  10'14 

Co 

1 

100 

58.93 

1.02  x  1019 

60Co 

100 

2772276  2.5  x  10'7 

15 

2.61  x  10'12 

Cu 

1 

69.17 

62.93 

9.57  x  1018 

64Cu 

100 

762.06 

0.00091 

15 

7.71  x  10'16 

Fe 

1 

0.282 

57.93 

1.04  x  1019 

59Fe 

100 

64080 

1.1  x  10-5 

15 

5.93  x  10'14 

*  Required  inputs 

Isotope  data  from  Chart  of  the  Nuclides,  16th  Ed. 

The  first  experiment  conducted  at  the  OSU  research  reactor  for  this  thesis  was  to 
characterize  the  flux  at  the  end  of  the  beam  port  with  the  collimation  apparatus  of 
polyethylene  plugs  and  cadmium  shielding  installed  as  shown  in  Figure  12.  The  outside 
of  the  beam  port  was  shielded  with  lead,  cadmium  and  boronated  bricks.  The  foils  were 
located  21.5  cm  outside  the  beam  port  as  shown  in  Figure  16.  The  reactor  was  brought 
up  to  20%  power  (100  kW)  and  radiation  surveys  conducted  to  ensure  regulatory 
compliance  and  safety.  The  escaping  neutron  and  gamma  radiation  in  the  reactor 
building  was  approximately  60  mrem/h,  high  enough  that  the  run  was  shut  down  after 
only  15  minutes.  For  this  experiment,  only  three  foils  were  used  as  inputs  to  the 
SAND-II  program,  both  gold  foils  and  the  bare  copper.  The  rest  did  not  have  measurable 
activity.  The  generated  spectrum  is  shown  in  Figure  17.  It  is  interesting  to  see  the 
impact  the  polyethylene  and  cadmium  have  when  comparing  this  spectrum  to  the  rabbit 
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Diff  Flux (n  cm-2  s-1  eV-1) 


tube  spectrum  in  Figure  9.  The  polyethylene  has  smoothed  the  spectrum  while  thermal 
neutron  attenuation  is  evident  from  the  cadmium  and  divergence  by  a  factor  of  eight. 


Figure  16.  Characterization  foils  at  end  of  beam  port. 


100%  Power 


Figure  17.  Neutron  energy  spectrum  at  end  of  OSU  collimated  beam  port. 
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Continued  operation  could  not  safely  be  maintained  at  a  high  enough  power  in 
this  configuration  to  produce  a  flux  capable  of  activating  the  foils  for  further 
experimentation,  to  include  a  run  using  N20.  Massive  shielding  had  been  applied  but  the 
escaping  radiation  was  still  too  large. 

In  order  to  minimize  escaping  radiation  while  maintaining  maximum  flux,  a 
polyethylene  plug  was  placed  in  the  outside  45  cm  of  the  beam  port  as  shown  in  Figure 
15.  The  plug  would  thennalize  neutrons  that  could  then  be  captured  by  the  cadmium  and 
boronated  bricks  shielding  the  outside  of  the  beam  port.  Lead  would  scatter  the  gammas. 

Characterization  foils  were  located  21.5  cm  from  the  beam  port.  Shielding  was 
built  up  from  the  floor  and  from  the  viewpoint  of  the  radiation  coming  out  of  the  beam 
port  consisted  of  a  cadmium  layer  followed  by  a  layer  of  lead  and  a  layer  of  boronated 
bricks  as  shown  in  Figure  18.  The  reactor  was  brought  up  to  10%  power  and  radiation 
surveys  begun.  The  gamma  levels  were  quite  low  but  the  neutron  levels  were  extremely 
high  (e.g.  600  rnrem/h)  all  around  the  beam  port  so  the  run  was  terminated  after  15 
minutes  of  operation.  Both  gold  characterization  foils  had  enough  activation  to  be 
counted  but  also  had  error  quoted  at  two  sigma  of  14-17%.  Their  calculated  saturated 
activities  per  nucleus  were  input  into  SAND-II  and  a  piecemeal  spectrum  came  out.  The 
integral  neutron  flux  was  approximately  10  n  cm'"  s'  ,  seven  times  the  flux  from  the 
simulation. 
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Figure  18.  Shielding  on  end  of  OSU  beam  port. 


The  results  of  the  previous  experiment  were  not  known  when  the  nitrous  bottle 
was  placed  on  the  end  of  the  beam  port  and  the  foils  placed  1.5  cm  from  the  cylinder  as 
shown  in  Figure  19.  The  shielding  was  supplemented  by  hanging  sheets  of  cadmium 
over  the  outside  of  the  shielding  to  absorb  thermal  neutrons.  The  reactor  was  brought  up 
to  only  5%  power  initially  so  safety  surveys  could  be  conducted.  Again,  the  neutron 
levels  were  unsafe  and  the  reactor  was  shut  down  after  only  a  few  minutes  of  operation. 
The  neutron  levels  were  approximately  600  mrem/h  near  the  containment  structure  all 
around  the  beam  port  but  declined  to  approximately  300  mrem/h  behind  the  shielding.  It 
was  postulated  that  the  polyethylene  plug  was  scattering  the  fast  and  epithennal  neutrons 
out  into  the  concrete  containment  structure  where  they  were  getting  further  scattered  out 
of  the  wall  and  into  the  room.  This,  in  effect,  made  the  entire  wall  a  source.  The 
shielding  was  effective  but  only  covered  the  beam  port  opening.  This  was  not  observed 
during  the  simulation  as  the  shielding  outside  the  beam  port  was  not  modeled  and  no 
detectors  were  placed  beyond  the  foil  location. 
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Figure  19.  N20  in  place  at  end  of  OSU  beam  port. 

The  gold  characterization  foils  were  the  only  ones  that  had  enough  activity  to  be 
counted.  The  foils  were  counted  for  approximately  39  hours  to  minimize  error  but  the 
activity  in  the  bare  foil  was  0.0472  disintegrations  per  second  with  25.15%  error  at  two 
sigma  and  the  cadmium-covered  foil  had  0.0304  disintegrations  per  second  with  44.29% 
error.  The  saturated  activity  per  nucleon  was  two  orders  of  magnitude  lower  than 
previous  activities  when  calculated  for  SAND-II  input.  SAND-II  did  converge  but  the 
resulting  spectrum  was,  not  surprisingly,  again  piecemeal.  The  flux  was 
1.7  x  1 0  n  cm"-  s"  ,  two  orders  of  magnitude  less  than  without  the  N20  and  in  the  same 
range  as  the  simulation.  These  results  have  too  much  error  associated  with  them  to  make 
definitive  comparisons  but  the  effect  of  the  N20  in  the  steel  cylinder  decreasing  flux  is 
readily  observed. 
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VI.  Conclusions  and  Recommendations 


Conclusions 

The  only  true  validation  of  the  methodology  at  ORNL  is  to  conduct  an  outside  test 
as  using  N20  to  replicate  large  volumes  of  air  did  not  give  results  to  send  to  ORNL.  An 
outside  test  allows  air  over  ground  transport  of  neutrons,  complete  and  measurable 
activation  of  foils,  ability  to  conduct  spectroscopy  of  the  foils,  and  an  independent  and 
valid  set  of  data  to  transfer  to  ORNL  for  validation. 

SCALE  can  be  very  complex.  Accuracy  can  be  improved  if  care  and  time  are 
taken  to  replicate  the  geometry  and  materials  involved  as  precisely  as  possible.  The 
initial  simulations  of  the  neutron  flux  at  the  outer  end  of  the  OSU  reactor  beam  port  were 
off  by  a  factor  of  nearly  one  thousand.  Simulating  the  core  as  a  cube  instead  of  a  point, 
using  more  accurate  geometry,  and  better  simulation  of  the  reactor  materials  brought  the 
simulated  neutron  flux  to  within  a  factor  of  seven  of  the  measured  neutron  flux.  More 
precise  simulation  of  the  actual  core  and  further  refinement  of  the  geometry  and  material 
should  produce  results  that  are  even  more  accurate. 

Only  a  small  portion  of  the  SCALE  program  was  used  in  this  thesis.  The  SAS4 
module  was  used  so  the  modeling  could  be  in  3-D.  It  is  not  very  easy  to  leam  but  the 
interactive  input  makes  it  easy  to  operate. 

Neutron  induced  reactions  in  gold  and  cobalt  that  were  not  originally  considered 
could  be  of  use.  The  (n,y)  thermal  reaction  of  gold  is  very  sensitive  and  has  a  2.69  day 
half-life.  Cobalt  has  an  (n,2n)  reaction  with  a  threshold  of  10.6  MeV  and  a  half-life  of  71 
days.  The  (n,p)  reaction  of  cobalt  has  a  threshold  of  0.8  MeV  and  a  half-life  of  44.5  days. 
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It  only  makes  sense  to  include  these  reactions  since  the  spectroscopy  will  be  done 
anyway. 

The  error  associated  with  the  activity  measurements  of  the  foils  must  be 
minimized.  The  use  of  wires  does  not  significantly  reduce  the  accuracy  as  the  difference 
in  solid  angle  between  a  point  and  a  wire  was  shown  to  be  negligible  for  the  spectroscopy 
system  used.  The  importance  of  a  good  calibration  was  driven  home  as  initial  results  had 
over  50%  error.  Time  and  resources  will  be  at  a  premium  immediately  following  a 
domestic  nuclear  event  and  it  was  shown  that  the  foils  could  be  measured  in  three  groups 
to  maximize  these  resources.  However,  caution  must  be  used  when  measuring  the  gold 
as  the  massive  activity  can  saturate  the  system  resulting  in  high  dead  time.  Some  foils, 
such  as  indium,  require  some  cooling  time  so  the  reactions  of  interest  can  be  more  readily 
detected. 

Using  N20  or  liquid  air  to  simulate  large  volumes  of  air  is  a  very  attractive 
proposition  however  the  concept  is  not  as  simple  as  first  thought.  A  possible  explanation 
is  that  the  angle  of  approach  between  neutrons  and  nucleons  is  different  at  different 
distances  from  the  source  [20].  This  would  change  the  probability  of  interaction  and  add 
an  unknown  factor  to  the  solution.  For  example,  point  A  in  air  in  Figure  20  has  a 
corresponding  point  B  in  N20.  Point  A  would  have  a  different  probability  of  interacting 
with  a  neutron  leaving  the  source  at  a  specific  angle  than  would  Point  B.  The  problem  of 
different  angles  of  interaction  at  depth  needs  to  be  resolved. 


47 


Source 


Foil 


Source 


B' 


N20 

1.53 

g/cc 


Foil 


Figure  20.  N2O  and  air  “equivalency”. 

The  OSU  research  reactor  has  many  capabilities.  Irradiating  experiments  on  the 
outside  of  the  beam  port  is  not  currently  one  of  them.  The  beam  port  points  directly 
toward  administrative  offices  and  limits  configuration  based  upon  safety.  Shielding  has 
to  be  moved  into  place  once  the  experiment  is  set  up  and  then  removed  to  access  the 
experiment.  This  can  be  very  time  consuming  and  also  increase  exposure  to  dangerous 
levels.  The  escaping  radiation  is  a  safety  concern,  especially  when  the  items  in  the  outer 
end  of  the  beam  port  scatter  neutrons  into  the  cement  containment  structure  and  then  into 
the  room.  A  substantial  shielding  plan  is  required  for  even  the  most  basic  operations. 


Recommendations 

The  use  of  the  OSU  reactor  beam  port  should  be  enhanced.  Easy-to-use  shielding 
should  be  designed  to  seal  around  the  open  beam  port  and  contain  escaping  radiation. 

One  idea  is  to  fabricate  a  tank  that  seals  around  the  port  leaving  a  cavity  for  experiments. 
The  tank  could  then  be  filled  with  ionized  water  that  is  available  at  OSU  to  provide  a 
moderator  for  the  neutrons.  Once  slowed  down,  the  neutrons  could  be  captured  by 
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cadmium  sheets.  Another  idea  is  to  build  shielding  on  a  pallet-type  jack.  The  shielding 
could  then  be  raised  and  lowered  hydraulically. 

I  would  not  recommend  pursuing  the  use  of  N20  or  liquid  air  at  OSU  unless  it 
could  be  placed  inside  the  beam  port.  Obtaining  approval  to  place  pressurized  containers 
in  the  beam  port  where  they  are  subject  to  heating  would  be  difficult.  Results  would  be 
better  if  the  foils  are  located  far  enough  away  from  the  N20  so  that  the  neutrons  coming 
through  the  N20  appear  as  from  a  point  source.  This  requires  space  and  shielding  along 
the  complete  path  between  the  N20  and  the  foil  to  keep  neutrons  from  scattering  in  from 
leakage  around  the  N20.  OSU  does  not  have  much  space  and  the  shielding  would  have 
to  contend  with  gammas  as  well  as  neutrons.  On  the  other  hand,  if  a  facility  with  a 
dedicated  experimentation  room  was  available,  pursuing  N20  or  liquid  air  use  could  be 
very  beneficial.  The  shielding  between  the  N20  and  the  foils  would  only  have  to  be  for 
thennal  neutrons  and  backscatter  onto  the  foil  would  be  minimized.  A  fabricated 
container  built  specifically  for  the  experimentation  area  could  be  sized  for  whatever  depth 
is  desired.  The  requirement  for  shielding  between  the  N20  and  foils  might  even  be  able 
to  be  done  away  with  if  the  container  was  designed  to  not  have  any  leakage  around  it. 

Testing  the  foils  with  an  outside  shot  to  get  results  to  send  to  ORNL  for 
methodology  verification  is  the  best  follow-on  research  for  this  project.  It  could  include 
more  options  for  wire  choice.  Gold,  vanadium  and  cobalt  are  all  rather  expensive.  Other 
wire  options  in  the  epithermal  range  to  better  cover  the  spectrum  could  also  be 
researched.  Design  and  testing  of  a  package  for  the  foils  could  be  researched,  especially 
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if  they  are  going  to  production.  If  outside  shots  become  available,  or  areas  with  more 
room  to  work  with,  research  on  using  N?0  or  liquid  air  could  be  beneficial. 
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Appendix  A.  Pu-Be  Source 


This  thesis  used  a  plutonium-beryllium  (Pu-Be)  source.  A  Pu-Be  source 
generates  the  majority  of  its  neutrons  from  (a,n)  reactions  with  some  contribution  from 
Be(n,2n)  and  Pu(n,f).  Another  small  contributor  is  the  inelastic  scattering  of  alpha 
particles  by  the  beryllium  that  yields  low  energy  neutrons  [2].  A  Pu-Be  source  initially 
has  an  increasing  neutron  yield  due  to  small  quantities  of  Pu  in  the  source.  Plutonium- 
241  is  a  beta  emitter  that  decays  to  24 ’Am.  Americium-241  is  an  alpha  emitter  which 
generates  more  (a,n)  reactions  resulting  in  more  neutrons.  After  several  decades,  the 
‘  Am  reaches  a  maximum  and  gradually  starts  to  decline. 

The  total  yield  of  a  stable  source  at  any  time  is  derived  by  M.  Edward  Anderson 
of  Mound  Laboratory  in  a  1968  Nuclear  Applications  article  [1].  It  is  given  as 


0(0  =  0(O){1  +  —  - —  («-'•'  -  )} 

a j  (TjAj  +  A2N ,  )(2,3  —  2,4) 


where 


Q(0)  =  alAlNl+ a2A2N2.  (16) 

The  time  prior  to  source  fabrication  when  the  americium  that  was  formed  during 
the  birth  of  the  plutonium  is  separated  from  the  plutonium  is  designated  as  t  =  0  .  The 
average  neutron  yield  per  alpha  particle  is  a  ,  A  is  the  disintegration  constant  and  N  is 
the  number  of  atoms  present  at  time  t  =  0 .  The  subscript  1  corresponds  to  Pu,  2  to 
240Pu,  3  to  241Pu  and  4  to  241  Am.  Since  the  alpha  particle  energies  of  239Pu  and  240Pu  are 
almost  identical,  ax  ~  a2 ,  and 


(9(0)  «  al(A1Ni  +  A2N2) . 


(17) 
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The  value  of  ax  is  not  known  but  the  yield  at  the  time  of  the  source  fabrication 


( 0(h)  )  is  typically  known  from  the  source  documentation.  Using  this  we  can  solve  for 


0(0)  = 


0(0 


i+*4 


A^A2N2 


(18) 


(e  Ai>'  -e  /hh) 


<7j  (A | N j  +  A2 N2  )(^3  —  A4 ) 


'Vl  ' 


and  then  solve  for  Q(t )  using  Equation  15.  The  value  of  —  is  1.23  ±  0.05  by 

a\ 

experimentation. 

The  Pu-Be  source  used  in  this  experiment,  M-l  170  produced  at  Mound 
Laboratory,  had  an  initial  neutron  emission  rate  ( 0(h)  )  of  9.04  x  106  n  s'1  on  9  March 
1962  [16].  It  was  composed  of  39.28  g  of  beryllium  and  76.36  g  of  plutonium.  The 
composition  of  the  plutonium  can  be  inferred  from  Anderson’s  article  to  be  93.4%  Pu, 
6.1%  240Pu  and  0.55%  241Pu  with  a  birth  date  ( t0 )  of  25  April  1961.  All  this  infonnation 

along  with  half-lives  was  used  to  create  the  spreadsheet  in  Table  10.  Equation  13  was 
implemented  to  calculate  the  number  of  atoms.  Once  the  information  is  set  up,  Equation 
18  and  Equation  15  are  solved  to  give  the  yield  as  of  the  date  input  by  the  user  ( t ).  To 
check  the  validity  of  the  calculations,  the  yield  for  22  October  1987  was  calculated.  It 
was  1.09  x  107  n  s'1,  within  1.9%  of  Whitworth’s  experimental  value  of  1.07  x  107  n  s'1  ± 
3%  for  the  same  source  [33]. 
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Table  10.  M-l  170  Pu-Be  source  calculator. 
Calculations  to  determine  activity  of  M-l  170  Pu-Be  source. 


N1  = 

1.79706  x  1023 

to  = 

4/25/1961 

N2  = 

1.16878  x  1022 

tl=(yr) 

3/9/1962 

0.871233 

N3  = 

1.04944  x  1021 

t  =  (yr) 

10/22/2004 

43.52329 

M  =  (/yr) 

2.84077  x  10'5 

X2  =  (/yr) 

0.000105022 

Q(tl)  =  (n/s) 

9.04  x  106 

k3  =  (/yr) 

0.05251115 

-U 

1! 

Tf- 

0.001513422 

Q(t0)  =  (n/s) 

8917066.151 

Q(0  =  (n/s) 

11280868.71 

Change  the  date  in  E5  to  get  the  Q(t)  of  that  date. 

Based  on  M.  Edward  Anderson's  article,  "Increases  in  Neutron  Yields  of  Plutonium-Beryllium  (a,n) 
Sources" 

Nuclear  Applications,  Volume  4,  March  1968 

The  Pu-Be  source  is  contained  in  a  barrel  of  paraffin  as  shown  in  Figure  7.  The 
paraffin  is  approximately  5  inches  thick  and  the  source  is  approximately  21.5  cm  from 
the  end  of  the  barrel.  The  yield  on  22  October  2004  when  irradiation  was  conducted  was 
calculated  to  be  1.128  x  107  n  s'1.  Using  the  theory  behind  Equation  1,  the  flux  at  the 
mouth  of  the  barrel  was  1.9  x  10  n  cm'"  s'  and  a  15-hour  irradiation  gave  a  fluence  on 
the  order  of  1  x  10  n  cm'  not  accounting  for  any  decay  during  the  activation. 

Neutron  energy  spectra  vary  among  different  Pu-Be  sources.  It  has  been  shown 
that  generally  larger  sources  have  a  higher  percentage  of  0-3  MeV  neutrons  [2].  Sources 
at  Mound  Laboratory  contained  between  1  and  160  g  of  plutonium.  Source  M-l  170 
contained  76.36  g.  A  spectrum  for  an  80  g  source  is  shown  in  Figure  21. 
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RELATIVE  NEUTRON  INTENSITY  (M eV*') 


1 


1 


NEUTRON  ENERGY  (MeV) 

Source:  Anderson,  Nuclear  Instalments  and  Methods,  1972 


Figure  21.  80  g  Pu-Be  source  neutron  energy  spectrum. 
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Appendix  B.  Gammas  from  Common  Neutron  Reactions 


Table  1 1  contains  a  list  of  common  neutron  initiated  reactions  for  the  foils 
selected  for  this  thesis.  It  lists  the  percentage  of  the  naturally  occurring  isotope,  mode  of 
decay,  percentage  of  decay,  the  energy  of  the  resulting  gamma  rays  and  their  associated 
branching  ratio  [17].  This  data  can  then  be  sorted  by  gamma  energy  to  screen  for 
gammas  that  are  similar  in  energy  and  might  cause  summing  problems  during  counting. 
A  feeling  of  the  magnitude  of  the  summing  problem  can  be  gained  by  looking  at  the  half- 
life,  percentage  of  natural  occurrence,  the  percentage  of  decay  and  the  branching 
ratio  [3]. 

The  indented  reactions  are  secondary  reactions  or  reactions  of  daughters.  The 
results  of  the  15-hour  Pu-Be  are  listed  and  whether  or  not  the  reaction  is  one  of  the 
AFIDS  reactions  of  interest. 
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Table  11.  Common  gammas  from  selected  foils. 

Reaction  T‘A  f''  .  Decay  %  Decay  '  branching  %  ROI 

natural  [keV]  with  Pu-Be 


27Al(n,y)28Al 

2.25m 

100 

p- 

27Al(n,2n)26Al 

7.1  x  105y 

100 

8 

27Al(n,2n)26Al 

7.1  x  105y 

100 

8 

27Al(n,p)27Mg 

9.45m 

100 

p- 

27Al(n,p)27Mg 

9.45m 

100 

p- 

27Al(n,a)24Na 

14.95h 

100 

p- 

27Al(n,a)24Na 

14.95h 

100 

p- 

197Au(n,y)198Au 

2.6952d 

100 

p- 

1 97  Au(n,y) 1 98  Aum 

2.27d 

100 

IT 

1 97  Au(n,y) 1 98  Aum 

2.27d 

100 

IT 

1 97  Au(n,y) 1 98  Aum 

2.27d 

100 

IT 

1 97  Au(n,y) 1 98  Aum 

2.27d 

100 

IT 

1 97  Au(n,y) 1 98  Aum 

2.27d 

100 

IT 

1 97  Au(n,2n)  196Au 

6.167d 

100 

p- 

197Au(n,2n)196Au 

6.167d 

100 

e 

1 97  Au(n,2n)  196Au 

6.167d 

100 

e 

1 97  Au(n,2n)  196Aum 

9.6h 

100 

IT 

1 97  Au(n,2n)  196Aum 

9.6h 

100 

IT 

1 97  Au(n,2n)  196Aum 

9.6h 

100 

IT 

1 97  Au(n,2n)  196Aum 

9.6h 

100 

IT 

1 97  Au(n,2n)  196Aum 

9.6h 

100 

IT 

1 97  Au(n,2n)  196Aum 

9.6h 

100 

IT 

197Au(n,p)I97Pt 

19.96h 

100 

P- 

197Au(n,p)197Pt 

19.96h 

100 

P- 

197Au(n,p)197Ptm 

1.59h 

100 

P- 

197Au(n,p)197Ptm 

1.59h 

100 

IT 

197Au(n,a)194Ir 

19. 3h 

100 

P- 

197Au(n,a)194Ir 

19. 3h 

100 

P- 

197Au(n,a)194Ir 

19. 3h 

100 

P- 

197Au(n,a)194Irm 

1 7  Id 

100 

IT 

197Au(n,a)194Irm 

1 7  Id 

100 

IT 

197Au(n,a)194Irm 

1 7  Id 

100 

P- 

197Au(n,a)194Irm 

1 7  Id 

100 

P- 

197Au(n,a)194Irm 

1 7  Id 

100 

P- 

197Au(n,a)194Irm 

1 7  Id 

100 

P- 

197Au(n,a)194Irm 

1 7  Id 

100 

P- 

100 

1778.85 

100 

N 

N 

100 

1129.67 

2.5 

N 

N 

100 

1808.65 

99.76 

N 

N 

100 

843.76 

71.8 

N 

N 

100 

1014.44 

28 

N 

N 

100 

1368.633 

100 

N 

Y 

100 

2754.028 

99.944 

N 

Y 

100 

411.802 

95.58 

Y 

N 

100 

97.21 

69.3 

N 

N 

100 

180.31 

50 

N 

N 

100 

204.1 

40.8 

N 

N 

100 

214.89 

77 

N 

N 

100 

333.82 

17 

N 

N 

7.2 

426.1 

6.6 

N 

Y 

92.8 

333.03 

22.9 

N 

Y 

92.8 

355.73 

87 

N 

Y 

100 

137.69 

1.3 

N 

Y 

100 

147.81 

43 

N 

Y 

100 

168.37 

7.7 

N 

Y 

100 

188.27 

37 

N 

Y 

100 

285.49 

4.3 

N 

Y 

100 

316.19 

2.9 

N 

Y 

100 

77.35 

17 

N 

N 

100 

191.437 

3.7 

N 

N 

3.3 

279 

2.4 

N 

N 

96.7 

346.5 

11.1 

N 

N 

100 

293.541 

2.5 

N 

N 

100 

328.448 

13.1 

N 

N 

100 

645.146 

1.18 

N 

N 

100 

84.288 

4.6 

N 

N 

100 

112.23 

8.3 

N 

N 

100 

111.7 

8.9 

N 

N 

100 

324 

2 

N 

N 

100 

328.5 

93 

N 

N 

100 

338.8 

55 

N 

N 

100 

390.8 

35 

N 

N 

56 


Reaction 

T!4 

% 

natural 

Decay 

%  Decay 

y  energy 
[keV] 

branching  % 

Visible 
with  Pu-Be 

ROI 

197Au(n,a)194Irm 

171d 

100 

P- 

100 

482.6 

97 

N 

N 

197Au(n,a)194Irm 

1 7  Id 

100 

P- 

100 

562.4 

70 

N 

N 

197Au(n,a)194Irm 

171d 

100 

P- 

100 

600.5 

62 

N 

N 

197Au(n,a)194Irm 

171d 

100 

P- 

100 

687.8 

59 

N 

N 

197Au(n,a)194Irm 

171d 

100 

P- 

100 

1011.8 

3.6 

N 

N 

59Co(n,y)60Co 

5.271y 

100 

P- 

100 

1173.228 

99.85 

N 

Y 

59Co(n,y)60Co 

5.271y 

100 

P- 

100 

1332.492 

99.9826 

N 

Y 

59Co(n,2n)58Co 

70.88d 

100 

8 

100 

810.7593 

99.45 

N 

N 

59Co(n,p)59Fe 

44. 5d 

100 

P- 

100 

142.651 

1.02 

N 

N 

59Co(n,p)59Fe 

44. 5d 

100 

P- 

100 

192.343 

3.08 

N 

N 

59Co(n,p)59Fe 

44. 5d 

100 

P- 

100 

1099.245 

56.5 

N 

N 

59Co(n,p)59Fe 

44. 5d 

100 

P- 

100 

1291.59 

43.2 

N 

N 

59Co(n,a)56Mn 

2.578h 

100 

P- 

100 

846.754 

98.9 

N 

N 

59Co(n,a)56Mn 

2.578h 

100 

P- 

100 

1810.72 

27.2 

N 

N 

59Co(n,a)56Mn 

2.578h 

100 

P- 

100 

2113.05 

14.3 

N 

N 

115In(n,n')115In 

4.485h 

95.71 

IT 

95 

336.241 

45.8 

N 

Y 

115In(n,y)116Inm 

54.2m 

95.71 

P- 

100 

138.326 

3.29 

N 

N 

115In(n,y)116Inm 

54.2m 

95.71 

P- 

100 

416.86 

27.7 

Y 

N 

115In(n,y)116Inm 

54.2m 

95.71 

P- 

100 

818.7 

11.5 

Y 

N 

115In(n,y)116Inm 

54.2m 

95.71 

P- 

100 

1097.3 

56.2 

Y 

N 

115In(n,y)116Inm 

54.2m 

95.71 

P- 

100 

1293.54 

84.4 

Y 

N 

115In(n,y)116Inm 

54.2m 

95.71 

P- 

100 

1507.4 

10 

Y 

N 

115In(n,y)116Inm 

54.2m 

95.71 

P- 

100 

1753.8 

2.46 

N 

N 

115In(n,y)116Inm 

54.2m 

95.71 

P- 

100 

2112.1 

15.5 

Y 

N 

115In(n,2n)114In 

1.198m 

95.71 

all 

N 

115In(n,p)115Cd 

2.228d 

95.71 

P- 

100 

260.896 

1.94 

N 

N 

115In(n,p)115Cd 

2.228d 

95.71 

P- 

100 

336.241 

45.9 

N 

N 

115In(n,p)115Cd 

2.228d 

95.71 

P- 

100 

492.351 

8.03 

N 

N 

115In(n,p)115Cd 

2.228d 

95.71 

P- 

100 

527.901 

27.5 

N 

N 

115In(n,a)112Ag 

3.13h 

95.71 

P- 

100 

606.7 

3.1 

N 

N 

115In(n,a)ll2Ag 

3.13h 

95.71 

P- 

100 

617.4 

43 

N 

N 

115In(n,a)112Ag 

3.13h 

95.71 

P- 

100 

692.7 

1.08 

N 

N 

115In(n,a)112Ag 

3.13h 

95.71 

P- 

100 

694.8 

3 

N 

N 

115In(n,a)ll2Ag 

3.13h 

95.71 

P- 

100 

851.2 

1.03 

N 

N 

ll5ln(n,a)ll2Ag 

3.13h 

95.71 

P- 

100 

1312.3 

1.2 

N 

N 

115In(n,a)112Ag 

3.13h 

95.71 

P- 

100 

1387.7 

5.4 

N 

N 

57 


Reaction 

T!4 

% 

natural 

Decay 

%  Decay 

y  energy 
[keV] 

branching  % 

Visible 
with  Pu-Be 

ROI 

ll3In(n,a)112Ag 

3.13h 

95.71 

P- 

100 

1613.6 

2.8 

N 

N 

115In(n,a)112Ag 

3.13h 

95.71 

P- 

100 

2106.2 

2.4 

N 

N 

115In(n,a)112Ag 

3. 13h 

95.71 

P- 

100 

2506.8 

1.08 

N 

N 

113In(n,y)114In 

1.198m 

4.29 

all 

N 

113In(n,2n)112In 

14.4m 

4.29 

e 

56 

606.4 

1.11 

N 

N 

113In(n,2n)112ln 

14.4m 

4.29 

e 

56 

617.1 

4.6 

N 

N 

113In(n,2n)112Inm 

20.8m 

4.29 

IT 

100 

156.4 

13.2 

N 

N 

113In(n,p)113Cd 

7.7  x  1015y 

4.29 

all 

N 

113In(n,a)110Ag 

24.6s 

4.29 

P- 

99.7 

657.5 

4.5 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

446.812 

3.62 

N 

N 

113In(n,a)ll0Agm 

249. 8d 

4.29 

P- 

98.64 

620.3553 

2.67 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

657.76 

94.3 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

677.6217 

10.56 

N 

N 

113In(n,a)ll0Agm 

249. 8d 

4.29 

P- 

98.64 

687.0091 

6.44 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

706.676 

16.33 

N 

N 

113In(n,a)ll0Agm 

249. 8d 

4.29 

P- 

98.64 

744.2755 

4.77 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

763.9424 

22.62 

N 

N 

113In(n,a)ll0Agm 

249. 8d 

4.29 

P- 

98.64 

818.0244 

7.34 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

884.6781 

72.7 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

937.485 

34.2 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

1384.2931 

24.9 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

1475.7792 

4.17 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

1505.028 

13.6 

N 

N 

113In(n,a)110Agm 

249. 8d 

4.29 

P- 

98.64 

1562.294 

1.244 

N 

N 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

328.762 

20.3 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

432.493 

2.9 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

487.021 

45.5 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

751.637 

4.33 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

815.772 

23.28 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

867.846 

5.5 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

919.55 

2.66 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

925.189 

6.9 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

1596.21 

95.4 

N 

Y 

139La(n,y)140La 

1.678d 

99.91 

P- 

100 

2521.4 

3.46 

N 

Y 

139La(n,2n)138La 

1.05  x  10ny 

99.91 

8 

65.6 

1435.795 

65.6 

N 

N 

139La(n,2n)138La 

1.05  x  10ny 

99.91 

P- 

34.4 

788.742 

34.4 

N 

N 

139La(n,p)139Ba 

1.396h 

99.91 

P- 

100 

165.8575 

24 

N 

N 
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Reaction 

T!4 

% 

natural 

Decay 

%  Decay 

y  energy 
[keV] 

branching  % 

Visible 
with  Pu-Be 

ROI 

139La(n,a)136Cs 

13.16(1 

99.91 

P- 

100 

66.881 

4.79 

N 

N 

139La(n,a)l36Cs 

13.16d 

99.91 

P- 

100 

86.36 

5.18 

N 

N 

139La(n,a)136Cs 

13.16d 

99.91 

P- 

100 

153.246 

5.75 

N 

N 

139La(n,a)l36Cs 

13.16d 

99.91 

P- 

100 

163.92 

3.39 

N 

N 

139La(n,a)136Cs 

13.16d 

99.91 

P- 

100 

176.602 

10 

N 

N 

139La(n,a)l36Cs 

13.16d 

99.91 

P- 

100 

273.646 

11.1 

N 

N 

139La(n,a)l36Cs 

13.16d 

99.91 

P- 

100 

340.547 

42.2 

N 

N 

139La(n,a)l36Cs 

13.16d 

99.91 

P- 

100 

818.514 

99.704 

N 

N 

139La(n,a)136Cs 

13.16d 

99.91 

P- 

100 

1048.073 

80 

N 

N 

139La(n,a)136Cs 

13.16d 

99.91 

P- 

100 

1235.362 

20 

N 

N 

138La(n,y)139La 

stable 

0.09 

N 

138La(n,2n)137La 

6  x  104y 

0.09 

no 

N 

138La(n,p)138Ba 

stable 

0.09 

N 

138La(n,a)135Csm 

53m 

0.09 

IT 

100 

787.2 

100 

N 

N 

138La(n,a)l35Csm 

53m 

0.09 

IT 

100 

846.1 

96 

N 

N 

55Mn(n,y)56Mn 

2.578h 

100 

P- 

100 

846.754 

98.9 

Y 

N 

55Mn(n,y)56Mn 

2.578h 

100 

P- 

100 

1810.72 

27.2 

Y 

N 

55Mn(n,y)56Mn 

2.578h 

100 

P- 

100 

2113.05 

14.3 

Y 

N 

55Mn(n,2n)54Mn 

312. Id 

100 

e 

100 

834.848 

99.976 

N 

Y 

55Mn(n,p)55Cr 

3.497m 

100 

all 

N 

55Mn(n,a)52V 

3.76m 

100 

P- 

100 

1434.06 

100 

N 

N 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

84.6808 

2.65 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

100.1065 

14.1 

N 

Y 

181Ta(n,y)l82Ta 

114.43d 

99.988 

P- 

100 

113.6725 

1.88 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

152.4308 

6.93 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

156.3876 

2.64 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

179.3945 

3.08 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

198.3532 

1.44 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

222.1096 

7.49 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

229.322 

3.63 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

264.0752 

3.61 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

1001.695 

2.07 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

1121.3008 

34.9 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

1189.0503 

16.2 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

1221.4066 

27 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

1231.0157 

11.44 

N 

Y 
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Reaction 

T'A 

% 

natural 

Decay 

%  Decay 

y  energy 
[keV] 

branching  % 

Visible 
with  Pu-Be 

ROI 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

1257.4185 

1.49 

N 

Y 

181Ta(n,y)182Ta 

114.43d 

99.988 

P- 

100 

1289.1561 

1.349 

N 

Y 

181Ta(n,y)182Tam 

15.8m 

99.988 

IT 

100 

146.785 

37.2 

N 

Y 

181Ta(n,y)182Tam 

15.8m 

99.988 

IT 

100 

171.586 

49 

N 

Y 

181Ta(n,y)182Tam 

15.8m 

99.988 

IT 

100 

184.951 

24.5 

N 

Y 

181Ta(n,y)182Tam 

15.8m 

99.988 

IT 

100 

318.4 

6.9 

N 

Y 

181Ta(n,2n)180Ta 

8.15h 

99.988 

e 

86 

93.4 

4.51 

Y 

N 

181Ta(n,p)181Hf 

42. 4d 

99.988 

P- 

100 

133.021 

43.3 

N 

N 

181Ta(n,p)181Hf 

42. 4d 

99.988 

P- 

100 

136.26 

5.85 

N 

N 

181Ta(n,p)181Hf 

42. 4d 

99.988 

P- 

100 

345.93 

15.12 

N 

N 

181Ta(n,p)181Hf 

42. 4d 

99.988 

P- 

100 

482.18 

80.5 

N 

N 

181Ta(n,a)178Lu 

28.5m 

99.988 

P- 

100 

93.179 

6 

N 

N 

181Ta(n,a)l78Lu 

28.5m 

99.988 

P- 

100 

1309.9 

1.4 

N 

N 

181Ta(n,a)178Lu 

28.5m 

99.988 

P- 

100 

1340.8 

3.42 

N 

N 

181Ta(n,a)178Lum 

23.1m 

99.988 

P- 

100 

88.85 

64.4 

N 

N 

181Ta(n,a)178Lum 

23.1m 

99.988 

P- 

100 

93.15 

17.2 

N 

N 

181Ta(n,a)178Lum 

23.1m 

99.988 

P- 

100 

213.41 

81.4 

N 

N 

181Ta(n,a)178Lum 

23.1m 

99.988 

P- 

100 

216.64 

2.48 

N 

N 

181Ta(n,a)178Lum 

23.1m 

99.988 

P- 

100 

325.6 

94.1 

N 

N 

181Ta(n,a)178Lum 

23.1m 

99.988 

P- 

100 

331.62 

11.4 

N 

N 

181Ta(n,a)178Lum 

23.1m 

99.988 

P- 

100 

426.36 

97 

N 

N 

51V(n,y)52V 

3.76m 

99.75 

P- 

100 

1434.06 

100 

N 

N 

51V(n,2n)50V 

1.4  x  1017y 

99.75 

P- 

17 

783.29 

17 

N 

N 

51V(n,2n)50V 

1.4  x  1017y 

99.75 

8 

83 

1553.77 

83 

N 

N 

51V(n,p)51Ti 

5.76m 

99.75 

P- 

100 

320.076 

93.1 

N 

N 

51V(n,p)51Ti 

5.76m 

99.75 

P- 

100 

608.55 

1.18 

N 

N 

51V(n,p)51Ti 

5.76m 

99.75 

P- 

100 

928.63 

6.9 

N 

N 

51V(n,a)48Sc 

43. 7h 

99.75 

P- 

100 

175.361 

7.48 

N 

Y 

51V(n,a)48Sc 

43. 7h 

99.75 

P- 

100 

983.526 

100.1 

N 

Y 

51V(n,a)48Sc 

43. 7h 

99.75 

P- 

100 

1037.522 

97.6 

N 

Y 

51V(n,a)48Sc 

43. 7h 

99.75 

P- 

100 

1212.88 

2.38 

N 

Y 

51V(n,a)48Sc 

43. 7h 

99.75 

P- 

100 

1312.12 

100.1 

N 

Y 

50V(n,y)51V 

stable 

0.25 

N 

50V(n,2n)49V 

33  Id 

0.25 

no 

N 

50V(n,p)50Ti 

stable 

0.25 

N 

50V(n,a)47Sc 

3.349d 

0.25 

p- 

100 

159.381 

68.3 

N 

N 
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Reaction 

T'A 

% 

natural 

Decay 

58Ni(n,y)59Ni 

7.6  x  104y 

68.08 

no 

58Ni(n,2n)57Ni 

35. 6h 

68.08 

p+ 

58Ni(n,2n)57Ni 

35. 6h 

68.08 

p+ 

58Ni(n,2n)57Ni 

35. 6h 

68.08 

p+ 

58Ni(n,2n)57Ni 

35. 6h 

68.08 

p+ 

57Ni  -  57Co  -  57Fe 

271. 8d 

- 

8 

57Ni  -  57Co  -  57Fe 

271. 8d 

- 

8 

58Ni(n,p)58Co 

70.88d 

68.08 

8 

58Ni(n,a)55Fe 

2.73y 

68.08 

no 

60Ni(n,y)61Ni 

stable 

26.22 

60Ni(n,2n)59Ni 

7.6  x  104y 

26.22 

no 

60Ni(n,p)60Co 

5.271y 

26.22 

p- 

60Ni(n,p)60Co 

5.271y 

26.22 

p- 

60Ni(n,a)57Fe 

stable 

26.22 

61Ni(n,y)62Ni 

stable 

1.14 

61Ni(n,2n)60Ni 

stable 

1.14 

61Ni(n,p)61Co 

1.65h 

1.14 

p- 

61Ni(n,p)61Co 

1.65h 

1.14 

p- 

61Ni(n,a)58Fe 

stable 

1.14 

62Ni(n,y)63Ni 

lOly 

3.63 

no 

62Ni(n,2n)61Ni 

stable 

3.63 

62Ni(n,p)62Co 

1.5m 

3.63 

p- 

62Ni(n,p)62Co 

1.5m 

3.63 

p- 

62Ni(n,p)62Co 

1.5m 

3.63 

p- 

62Ni(n,p)62Co 

1.5m 

3.63 

p- 

62Ni(n,p)62Co 

1.5m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,p)62Co 

13.9m 

3.63 

p- 

62Ni(n,a)59Fe 

44. 5d 

3.63 

p- 

62Ni(n,a)59Fe 

44. 5d 

3.63 

p- 

Decay 

y  energy 
[keV] 

branching  % 

Visible 
with  Pu-Be 

ROI 

N 

100 

127.164 

16.7 

N 

Y 

100 

1377.63 

81.7 

N 

Y 

100 

1757.55 

5.75 

N 

Y 

100 

1919.52 

12.3 

N 

Y 

100 

122.06065 

85.6 

N 

Y 

100 

136.47356 

10.68 

N 

Y 

100 

810.7593 

99.45 

N 

Y 

N 

N 

N 

100 

1173.228 

99.85 

N 

N 

100 

1332.492 

99.9826 

N 

N 

N 

N 

N 

100 

67.415 

84.7 

N 

N 

100 

909.2 

3.6 

N 

N 

N 

N 

N 

100 

1128.9 

11.1 

N 

N 

100 

1172.9 

83.5 

N 

N 

100 

1985 

1.6 

N 

N 

100 

2301.8 

14.7 

N 

N 

100 

2345.9 

1.3 

N 

N 

99 

777.5 

1.72 

N 

N 

99 

875 

1.28 

N 

N 

99 

1129 

1.29 

N 

N 

99 

1163.5 

68 

N 

N 

99 

1172.9 

97.7 

N 

N 

99 

1718.7 

6.7 

N 

N 

99 

2003.7 

18.6 

N 

N 

99 

2104.7 

6.4 

N 

N 

99 

2301.9 

1.79 

N 

N 

99 

2882.3 

1.09 

N 

N 

100 

142.651 

1.02 

N 

N 

100 

192.343 

3.08 

N 

N 
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Reaction 

T!4 

% 

natural 

Decay 

%  Decay 

y  energy 
[keV] 

branching  % 

Visible 
with  Pu-Be 

ROI 

62Ni(n,a)59Fe 

44. 5d 

3.63 

P- 

100 

1099.245 

56.5 

N 

N 

62Ni(n,a)59Fe 

44. 5d 

3.63 

P- 

100 

1291.59 

43.2 

N 

N 

64Ni(n,y)65Ni 

2.517h 

0.93 

P- 

100 

366.27 

4.81 

N 

N 

64Ni(n,y)65Ni 

2.517h 

0.93 

P- 

100 

1115.53 

15.43 

N 

N 

64Ni(n,y)65Ni 

2.517h 

0.93 

P- 

100 

1481.84 

23.59 

N 

N 

64Ni(n,2n)63Ni 

lOly 

0.93 

no 

N 

64Ni(n,p)64Co 

0.3s 

0.93 

p- 

100 

931.1 

5 

N 

N 

64Ni(n,p)64Co 

0.3s 

0.93 

p- 

100 

1346.1 

10 

N 

N 

64Ni(n,a)6lFe 

6m 

0.93 

p- 

100 

120.34 

5.3 

N 

N 

64Ni(n,a)6lFe 

6m 

0.93 

p- 

100 

177.61 

2 

N 

N 

64Ni(n,a)61Fe 

6m 

0.93 

p- 

100 

297.9 

22 

N 

N 

64Ni(n,a)61Fe 

6m 

0.93 

p- 

100 

1027.42 

43 

N 

N 

64Ni(n,a)61Fe 

6m 

0.93 

p- 

100 

1205.07 

44 

N 

N 

64Ni(n,a)6lFe 

6m 

0.93 

p- 

100 

1645.95 

7 

N 

N 

64Ni(n,a)61Fe 

6m 

0.93 

p- 

100 

2011.6 

4.4 

N 

N 

61Fe-61Co-6lNi 

1.65h 

1.14 

p- 

100 

67.415 

84.7 

N 

N 

61Fe  -  61Co  -  6lNi 

1.65h 

1.14 

p- 

100 

909.2 

3.6 

N 

N 

63Cu(n,y)64Cu 

12.70  lh 

69.17 

all 

Y 

63Cu(n,2n)62Cu 

9.74m 

69.17 

all 

N 

63Cu(n,p)63Ni 

lOly 

69.17 

no 

N 

63Cu(n,a)60Co 

5.271y 

69.17 

P- 

100 

1173.228 

99.85 

N 

Y 

63Cu(n,a)60Co 

5.271y 

69.17 

P- 

100 

1332.492 

99.9826 

N 

Y 

65Cu(n,y)66Cu 

5.1m 

30.83 

P- 

100 

1039.2 

9.23 

N 

N 

65Cu(n,2n)64Cu 

12.70  lh 

30.83 

all 

Y 

65Cu(n,p)65Ni 

2.517h 

30.83 

P- 

100 

366.27 

4.81 

N 

N 

65Cu(n,p)65Ni 

2.517h 

30.83 

P- 

100 

1115.53 

15.43 

N 

N 

65Cu(n,p)65Ni 

2.517h 

30.83 

P- 

100 

1481.84 

23.59 

N 

N 

65Cu(n,a)62Co 

1.5m 

30.83 

P- 

100 

1128.9 

11.1 

N 

N 

65Cu(n,a)62Co 

1.5m 

30.83 

P- 

100 

1172.9 

83.5 

N 

N 

65Cu(n,a)62Co 

1.5m 

30.83 

P- 

100 

1985 

1.6 

N 

N 

65Cu(n,a)62Co 

1.5m 

30.83 

P- 

100 

2301.8 

14.7 

N 

N 

65Cu(n,a)62Co 

1.5m 

30.83 

P- 

100 

2345.9 

1.3 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

777.5 

1.72 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

875 

1.28 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

1129 

1.29 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

1163.5 

68 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

1172.9 

97.7 

N 

N 

62 


Reaction 

TVi 

% 

natural 

Decay 

%  Decay 

y  energy 
[keV] 

branching  % 

Visible 
with  Pu-Be 

ROI 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

1718.7 

6.7 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

2003.7 

18.6 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

2104.7 

6.4 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

2301.9 

1.79 

N 

N 

65Cu(n,a)62Co 

13.9m 

30.83 

P- 

99 

2882.3 

1.09 

N 

N 
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Appendix  C.  Point-to-Wire  Solid  Angle  Comparison 


The  accurate  counting  of  gamma  rays  has  a  major  dependence  on  the  geometry  of 
the  counting  setup.  The  absolute  efficiency  of  the  counting  system  is  the  ratio  of  the 
number  of  gamma  rays  counted  to  the  number  of  gamma  rays  emitted  by  the  source.  One 
factor  of  absolute  efficiency  is  the  solid  angle  subtended  by  the  detector  at  the  source 
position. 

For  a  point  source  located  along  the  axis  of  a  right  circular  cylindrical  detector, 
Figure  22,  the  solid  angle  (in  steradians)  is  given  by 

n  =  2*(i --=£=)  (19) 

yld2  +  a2 

where  d  is  the  source-detector  distance  and  a  is  the  detector  radius  [12]. 


Figure  22.  Solid  angle  for  a  point  source. 

The  HPGe  Canberra  Model  GC 10021  detector  has  a  crystal  that  is  83  mm  in 
diameter  so  a  is  41.5  mm.  Plotting  Equation  19  while  varying  d  from  0  to  100  mm 
yields  the  curve  in  Figure  23.  The  plot  passes  the  common  sense  test  as  it  starts  large  and 
tapers  off  as  the  distance  increases. 
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Figure  23.  Point  source  solid  angle  vs.  distance  based  on  Knoll  [12]. 


Source 


Detector 


Figure  24.  Solid  angle  geometry  for  a  circular  disk. 

This  thesis  used  wires  instead  of  point  sources.  No  equations  could  be  found  for 
wire  geometry  but  several  could  be  found  for  the  more  common  disk  geometry  shown  in 
Figure  24.  A  wire  shows  the  detector  a  smaller  surface  than  a  disk  with  a  diameter  the 
same  dimension  as  the  length  of  the  wire  so  the  solid  angle  must  also  be  smaller.  Using 
the  disk  as  a  worst-case  scenario  and  s  as  the  radius  of  the  source,  it  has  been  shown  that 


2 m  r  e  dkJ j  ( sk)Jx  ( ak ) 


Q=^f 

c  J 
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where  the  Jx(x)  are  Bessel  functions  of  x  [12].  Since  Equation  20  does  not  have  an 


analytic  solution,  Knoll  gives  an  approximation  as  shown  in  Figure  25  where  the 
variables  are  as  shown  in  Figure  24. 

=.  (*)2 

S\2 


@[d_]  (-) 


Fi[d  ]  : 


d ' 
5 


fii  d] 


35 


/3  [d]  ■ 


F2[d  ]  : 


-  16  (1  +  /3  [d] ) 7/2  16  (1  +  /3  [d] ) 9/2 

35  /3[d]  315  /3  [d] 2 
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/3[d] 


-  128  (1  +  /3  [d] )  9/2  256  (1  +  /3  [d] ) u/2  1024  (1  +  /3  [d] ) 13/2 

*fl-  1 


nc[d  ]  :=  2 


—  -  —  a[d]  ^[d]  +  a[d]2  Fi  [d]  -a[d]3F2[d]l 

(l  +  /3[d])1/2  8  (1  +  /3  [d] ) 5/2  ; 


Figure  25.  Approximation  of  solid  angle  for  disk  geometry  based  on  Knoll  [12]. 
Since  the  majority  of  wires  used  in  this  thesis  are  1  cm  long,  the  radius  of  the  disk 
source  is  set  at  5  mm.  The  radius  of  the  detector  remains  at  41 .5  mm.  The  plot  of  the 
approximation  equation  in  Figure  25  for  distances  from  40  to  100  mm  is  given  in 
Figure  26.  It  is  obvious  that  this  approximation  is  not  accurate  at  shorter  distances  but 


appears  to  closely  resemble  the  solid  angle  curve  in  Figure  23  from  about  55  mm  to  at 
least  100  mm.  For  the  detector  used  in  this  thesis,  the  crystal  is  5  mm  from  the  face  of 
the  detector  and  the  wires  are  measured  1  mm  above  that  for  a  total  of  6  mm.  The 


approximation  does  not  work  for  this  situation. 
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Figure  26.  Approximation  of  disk  solid  angle  vs.  distance  based  on  Knoll  [12]. 
Ruby  gives  another  approximation  for  disk  geometry  as  shown  in  Figure  27  [21]. 
The  variables  are  the  same  as  in  Knoll.  Again,  it  is  not  accurate  near  the  face  of  the 
detector  but  is  a  close  approximation  from  about  40  mm  to  at  least  100  mm  as  shown  in 
Figure  28.  Again,  this  does  not  work  with  the  distance  of  6  mm. 

«[dJ  :■  (|)2 

.[dJ  :■  (|)2 

Qx[d_]  := 

3  15 

7r*a[d]  (l--  03 [d]  -a[d])  +  —  (/3[d]2  +  a[d]2  +  3j8[d]  a[ d] )  - 
35 

—  (£[d]3  +  a[d]3+  6/3  [d]  a[d]  (J8[d]  +  at [d] ) ) ) 

Figure  27.  Approximation  of  solid  angle  for  disk  geometry  based  on  Ruby  [21]. 


67 


Figure  28.  Approximation  of  disk  solid  angle  vs.  distance  based  on  Ruby  [21]. 

The  integral  of  the  Bessel  functions  in  Equation  20  might  not  have  analytic 
solutions  but  Mathematica  can  easily  approximate  it  using  numerical  techniques. 
Values  were  calculated  using  Mathematica  and  Equations  19  and  20  from  1  mm  to 
1000  mm  and  the  results  tabulated  in  Table  12.  Again,  the  detector  radius  was  41.5  mm 
and  the  source  radius  was  5  mm.  The  error  maximized  at  about  35  mm.  Closer  than 
35  mm  the  size  of  the  detector  overpowers  the  1  cm  disk  and  over  35  mm  the  disk 
appears  increasingly  like  a  point  to  the  detector.  At  1  mm  off  the  surface  of  the  detector 
where  the  wires  are  measured  (6  mm  overall),  the  difference  in  solid  angle  between  a 
point  and  a  1  cm  disk  is  less  than  0.1%.  Since  a  wire  has  a  smaller  solid  angle  than  a 
disk,  the  difference  would  be  even  less.  This  amount  of  error  is  negligible  when 
compared  to  the  statistical  error  of  spectroscopy. 
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Table  12.  Point  vs.  disk  solid  angle  comparison. 

Comparison  of  Knoll  (pg  118)  point  and  disk  solid  angle  equations 
Calculations  done  in  Mathematica 

Detector  radius  =  41.5  mm  Source  radius  =  5  mm 


Distance  (mm) 

Solid  Angle 

disk 

point 

%  difference 

1 

6.1310092 

6.1318271 

0.01% 

5 

5.5276006 

5.5316102 

0.07% 

6 

5.3793851 

5.3841208 

0.09% 

7 

5.2327099 

5.2381331 

0.10% 

8 

5.0877981 

5.093865 

0.12% 

9 

4.944858 

4.951522 

0.13% 

10 

4.8040826 

4.8112934 

0.15% 

15 

4.1382386 

4.1473868 

0.22% 

20 

3.5455719 

3.5553911 

0.28% 

25 

3.0314569 

3.0409823 

0.31% 

30 

2.5935465 

2.6022 

0.33% 

35 

2.2248674 

2.2323981 

0.34% 

40 

1.9164368 

1.9228116 

0.33% 

50 

1.4440195 

1.4483921 

0.30% 

60 

1.1127091 

1.1156496 

0.26% 

70 

0.8764548 

0.8784413 

0.23% 

80 

0.7044239 

0.7057867 

0.19% 

90 

0.5764273 

0.5773808 

0.17% 

100 

0.4792129 

0.4798939 

0.14% 

102 

0.4626302 

0.4632684 

0.14% 

125 

0.3197317 

0.3200513 

0.10% 

150 

0.2273252 

0.2274918 

0.07% 

200 

0.1309911 

0.1310482 

0.04% 

250 

0.0847964 

0.0848207 

0.03% 

500 

0.0215296 

0.0215312 

0.01% 

1000 

0.0054035 

0.0054036 

0.00% 
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Appendix  D.  Wire  Activation  Simulation  Results 


Table  13.  Calculations  to  determine  wire  size. 


(cj(|>)  SCALE 

density 

radius 

length 

volume 

mass 

#  of  atoms 

#  of  atoms 

activity 

T  !4 

X 

A  at  +4  days 

Al 

500  m 

[n/s] 

[g/crrv3] 

[cm] 

[cm] 

[cm3] 

[mg] 

target 

activated 

[Bq] 

[days] 

[/d] 

[Bq] 

6.26x1  O'13 

2.70 

0.05 

1.00 

7.85  x  1 0'3 

21.20 

4.73  x  102“ 

2.96  x  10“ 

3815.97 

0.62 

1.11 

4.46x10’ 

330  m 

3.24x1 0'1Z 

2.70 

0.05 

0.05 

3.71  x  10" 

1.00 

2.23  x  1019 

7.24  xIO7 

932.27 

0.62 

1.11 

1.09x10’ 

700  m 

2.07x1  O'13 

2.70 

0.05 

1.00 

7.85  x  1 0'3 

21.20 

4.73  x  102“ 

9.78  x  107 

1259.02 

0.62 

1.11 

1.47x10’ 

800  m 

7.64x10'" 

2.70 

0.05 

1.00 

7.85  x  1 0'3 

21.20 

4.73  x  102“ 

3.61  x  107 

465.41 

0.62 

1.11 

5.44 

1000  m 

1.84x10'" 

2.70 

0.05 

1.00 

7.85  x  1 0'3 

21.20 

4.73  x  102“ 

8.69  x  10“ 

111.95 

0.62 

1.11 

1.31 

2500  m 

Au 

500  m 

5.78  x  10’’“ 

2.70 

0.05 

1.00 

7.85  x  1 0'3 

21.20 

4.73  x  102“ 

2.73  x  102 

0.00 

0.62 

1.11 

4.11  x  10'“ 

9.46x1  O'12 

19.30 

0.05 

0.30 

2.36  x  1 0'3 

45.47 

1 .39  x  102“ 

1.32  x  10“ 

1707.56 

6.18 

1.12x10" 

1.09  x  10“ 

1000  m 

2.52x1  O'13 

19.30 

0.05 

0.30 

2.36  x  1 0'3 

45.47 

1 .39  x  102“ 

3.50  x  107 

45.47 

6.18 

1.12x10" 

2.90  x  10’ 

1100  m 

1.85x1  O'13 

19.30 

0.05 

0.30 

2.36  x  1 0'3 

45.47 

1 .39  x  102“ 

2.57  x  107 

33.33 

6.18 

1.12x10" 

2.13x  10’ 

1200  m 

6.22x10" 

19.30 

0.05 

0.30 

2.36  x  1 0'3 

45.47 

1 .39  x  102“ 

8.65  x  10“ 

11.22 

6.18 

1.12x10" 

7.17 

2500  m 

Co 

500  m 

4.91  x  1 0'18 

19.30 

0.05 

0.30 

2.36  x  1 0'3 

45.47 

1 .39  x  102“ 

6.82  x  102 

0.00 

6.18 

1.12x  10" 

5.65  x  10" 

6.24x10'" 

8.90 

0.05 

1.00 

7.85  x  1 0'3 

69.90 

7.14  x  102“ 

4.46x10™ 

185.90 

1923.92 

3.60x10" 

1.86x10“ 

1000  m 

1.05x10'" 

8.90 

0.05 

1.00 

7.85  x  1 0'3 

69.90 

7.14  x  102“ 

7.50  x  10“ 

31.25 

1923.92 

3.60x10" 

3.12x  10’ 

1100  m 

3.80x1  O'12 

8.90 

0.05 

1.00 

7.85  x  1 0'3 

69.90 

7.14  x  102“ 

2.72  x  10“ 

11.33 

1923.92 

3.60x10" 

1.13x  10’ 

1200  m 

8.38x1  O'13 

8.90 

0.05 

1.00 

7.85  x  1 0'3 

69.90 

7.14  x  102“ 

5.99  x  10“ 

2.50 

1923.92 

3.60x10" 

2.49 

2500  m 

In 

100  m 

1.30x1  O'15 

8.90 

0.05 

1.00 

7.85  x  1 0'3 

69.90 

7.14  x  102“ 

9.27  x  10“ 

0.00 

1923.92 

3.60x10" 

3.86  x  1 0'3 

1.22  x  10  “ 

7.31 

0.05 

1.00 

7.85  x  1 0'3 

57.41 

3.01  x  102“ 

3.69x10" 

1.58  xIO7 

0.19 

3.71 

5.73 

200  m 

2.89x10"“ 

7.31 

0.05 

1.00 

7.85  x  1 0'3 

57.41 

3.01  x  102“ 

8.70x10™ 

3.73  xIO6 

0.19 

3.71 

1.35 

300  m 

1.09x10"“ 

7.31 

0.05 

1.00 

7.85  x  1 0'3 

57.41 

3.01  x  102“ 

3.28x10™ 

1.41  xIO6 

0.19 

3.71 

5.10  x  10" 

500  m 

2.49x10'" 

7.31 

0.05 

1.00 

7.85  x  1 0'3 

57.41 

3.01  x  102“ 

7.50  x  10“ 

3.22  xIO5 

0.19 

3.71 

1.17x  10" 

1000  m 

1.75  xIO"2 

7.31 

0.05 

1.00 

7.85  x  1 0'3 

57.41 

3.01  x  102“ 

5.27  x  10“ 

2.26  xIO4 

0.19 

3.71 

8.20  x  1 0’3 

2500  m 

La 

500  m 

7.62  xIO"7 

7.31 

0.05 

1.00 

7.85  x  1 0'3 

57.41 

3.01  x  102“ 

2.29  x  104 

0.98 

0.19 

3.71 

3.57  x  10  7 

1.15x10'" 

6.15 

0.05 

0.50 

3.93  x  1 0'3 

24.13 

1 .05  x  102“ 

1.20  x  10“ 

5747.28 

1.68 

4.13x10" 

1 . 1 0  x  10“ 

1000  m 

5.31  xIO"2 

6.15 

0.05 

0.50 

3.93  x  1 0'3 

24.13 

1 .05  x  102“ 

5.55  x  10“ 

2654.46 

1.68 

4.13x10" 

5.09  x  10“ 

1600  m 

2.72x10"“ 

6.15 

0.05 

0.50 

3.93  x  1 0'3 

24.13 

1 .05  x  102“ 

2.85  x  107 

136.26 

1.68 

4.13x10" 

2.61  x  10’ 

1700  m 

1.32x10'" 

6.15 

0.05 

0.50 

3.93  x  1 0'3 

24.13 

1 .05  x  102“ 

1.39  x  10“ 

6.62 

1.68 

4.13  x  10" 

1.27 

2500  m 

Mn  80% 

1.29x10'" 

6.15 

0.05 

0.50 

3.93  x  1 0'3 

24.13 

1 .05  x  102“ 

1.35  x  10“ 

6.48 

1.68 

4.13  x  10" 

1.24 

500  m 

2.71  xIO"2 

7.30 

0.05 

1.50 

1.18  x  10  '2 

68.80 

7.54  x  102“ 

2.05  x  10“ 

52.60 

312.10 

2.22x10'“ 

5.21  x  10’ 

700  m 

5.23x10"“ 

7.30 

0.05 

1.50 

1.18  x  10  '2 

68.80 

7.54  x  102“ 

3.95  x  10“ 

10.14 

312.10 

2.22x10'“ 

1.01  x  10’ 

1000  m 

4.56x10'" 

7.30 

0.05 

1.50 

1.18  x  10  '2 

68.80 

7.54  x  102“ 

3.44  xIO7 

0.88 

312.10 

2.22x10'“ 

8.75  x  10" 

1200  m 

7.90x1  O'15 

7.30 

0.05 

1.50 

1.18  x  10z 

68.80 

7.54  x  102“ 

5.96  x  10“ 

0.15 

312.10 

2.22x10'“ 

1.52  x  10" 

2500  m 

Ta 

500  m 

3.61  x  1 0'18 

7.30 

0.05 

1.50 

1.18  x  10  '2 

68.80 

7.54  x  102“ 

2.73  x  103 

0.00 

312.10 

2.22x10'“ 

6.94x10“ 

3.83x10"“ 

16.65 

0.02 

0.50 

5.67  x  10" 

9.44 

3.14  x  10'“ 

1.21  xIO™ 

844.82 

114.43 

6.06x10'“ 

8.25  x  10“ 

1000  m 

6.65x10'" 

16.65 

0.02 

0.50 

5.67  x  10" 

9.44 

3.14  x  10'“ 

2.09  x  10“ 

146.42 

114.43 

6.06x10'“ 

1.43  x  10“ 

1100  m 

1.33x10'" 

16.65 

0.02 

0.50 

5.67  x  10" 

9.44 

3.14  x  10'“ 

4.18  x  10“ 

29.33 

114.43 

6.06x10'“ 

2.86  x  10’ 

1200  m 

7.73  xIO"2 

16.65 

0.02 

0.50 

5.67  x  10" 

9.44 

3.14  x  10'“ 

2.43  x  10“ 

17.04 

114.43 

6.06x10'“ 

1.66x  10’ 

1300  m 

5.36  xIO"2 

16.65 

0.02 

0.50 

5.67  x  10" 

9.44 

3.14  x  10'“ 

1.68  x  10“ 

11.80 

114.43 

6.06x10'“ 

1.15x  10’ 

2500  m 

6.82x10"“ 

16.65 

0.02 

0.50 

5.67  x  10" 

9.44 

3.14  x  10'“ 

2.14x10’ 

0.00 

114.43 

6.06x10'“ 

1.47  x  10  “ 

V 

500  m 

6.23x10" 

6.11 

0.05 

1.00 

7.85  x  1 0'3 

47.99 

5.67  x  102“ 

3.53  x  107 

155.61 

1.82 

3.81  xIO" 

3.40  x  10’ 

600  m 

3.38x10" 

6.11 

0.05 

1.00 

7.85  x  1 0'3 

47.99 

5.67  x  102“ 

1.92  x  107 

84.41 

1.82 

3.81  xIO" 

1.84x  10’ 

700  m 

1.60x10'" 

6.11 

0.05 

1.00 

7.85  x  1 0'3 

47.99 

5.67  x  102“ 

9.08  x  10“ 

40.00 

1.82 

3.81  xIO" 

8.73 

1000  m 

1.75x10"“ 

6.11 

0.05 

1.00 

7.85  x  1 0'3 

47.99 

5.67  x  102“ 

9.95  x  10“ 

4.38 

1.82 

3.81  x  10" 

9.57  x  10" 

2500  m 

1 .58  x  1 0’2“ 

6.11 

0.05 

1.00 

7.85  x  1 0'3 

47.99 

5.67  x  102“ 

8.96 

0.00 

1.82 

3.81  x  10" 

8.61  x  10" 
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(ct<|>)  SCALE 

density 

radius 

length 

volume 

mass 

#  of  atoms 

#  of  atoms 

activity 

T  !4 

X 

A  at  +4  days 

Ni 

[n/s]  [g/cm3] 

cross-section  for  (n,p) 

[cm] 

[cm] 

[cm3] 

[mg] 

target 

activated 

[Bq] 

[days] 

[/d] 

[Bq] 

500  m 

3.67x1 0'1Z 

8.90 

0.05 

1.00 

7.85  x  icr3 

69.92 

7.17  x  10™ 

2.63  x  10s 

298.12 

70.88 

9.78  x  1 0'3 

2.87  x  102 

1000  m 

1.57x1  O'13 

8.90 

0.05 

1.00 

7.85  x  icr3 

69.92 

7.17  x  10™ 

1.13x  10“ 

12.74 

70.88 

9.78  x  1 0'3 

1.22  x  101 

1100  m 

9.82x10'" 

8.90 

0.05 

1.00 

7.85  x  icr3 

69.92 

7.17  x  10™ 

7.05  x  107 

7.97 

70.88 

9.78  x  1 0'3 

7.67 

2500  m 

Cu 

7.22x1 0'1S  8.90 

cross-section  for  (n,a) 

0.05 

1.00 

7.85  x  icr3 

69.92 

7.17  x  10™ 

5.18  x  103 

0.00 

70.88 

9.78  x  1 0'3 

5.63  x  10" 

250  m 

1.61  xIO'12 

8.96 

0.05 

2.00 

1.57  x  icr2 

140.74 

1.33  x  1021 

2.14  x10s 

8.93 

1923.92 

3.60  x  1 0'4 

8.92 

300  m 

9.69x1  O'13 

8.96 

0.05 

2.00 

1.57  x  10'2 

140.74 

1.33  x  1021 

1.29  x  10s 

5.39 

1923.92 

3.60  x  1 0'4 

5.38 

500  m 

1.51  xIO'13 

8.96 

0.05 

2.00 

1.57  x  icr2 

140.74 

1.33  x  1021 

2.02  x  10“ 

0.84 

1923.92 

3.60  x  1 0'4 

8.40x10" 

1000  m 

2.70x1  O'15 

8.96 

0.05 

2.00 

1.57  x  10'2 

140.74 

1.33  x  1021 

3.60  x  10“ 

0.02 

1923.92 

3.60  x  1 0'4 

1.50  x  1 0'2 

2500  m 

8.70x1  O'20 

8.96 

0.05 

2.00 

1.57  x  icr2 

140.74 

1.33  x  1021 

1.16  x  102 

0.00 

1923.92 

3.60  x  1 0'4 

4.83  x  10" 

Table  14.  MCNP  calculations  to  verify  SCALE  at  500  meters. 


Energy 

MCNP 

Initial 

Total  Flux 

Cross-section 

Convert  to 

Flux  times 

Target  Atoms 

atoms 

Group  [MeV] 

[/cm2] 

Neutrons 

[n/cm2] 

SCALE  [b] 

[cm2] 

cross-section 

lmg 

activated 

14 

1.758  x  1  O'12 

1.000  x 

1024 

1.758  x  1012 

1.226  x  104 

1.223  x  10'25 

2.155  x  1043 

2.233  x  1019 

4.812  x  10' 

13.84 

1.938  x  10'12 

1.000  x 

1024 

1.938  x  1012 

1.279  x  104 

1.279  x  10'25 

2.478  x  1043 

2.233  x  1019 

5.533  x  10' 

12.84 

2.023  x  1042 

1.000  x 

1024 

2.023  x  1012 

1.1 16  x  104 

1.1 16  x  10'25 

2.257  x  1043 

2.233  x  1019 

5.040  x  10' 

10 

6.342  x  1043 

1.000  x 

1024 

6.342  x  10" 

6.643  x  10'2 

6.643  x  10'26 

4.213  x  1044 

2.233  x  1019 

9.407  x  10' 

8.1873 

9.948  x  1043 

1.000  x 

1024 

9.948  x  10“ 

1.976  x  10'2 

1.976  x  10'26 

1.966  x  1044 

2.233  x  1019 

4.389  x  10' 

6.434 

1.734  x  1042 

1.000  x 

1024 

1.734  x  1012 

7.053  x  10*4 

7.053  x  10'28 

1.223  x  1045 

2.233  x  1019 

2.732  x  10‘ 

4.8 

6.904  x  1043 

1.000  x 

1024 

6.904  x  10" 

1.521  x  10'7 

1.521  x  10‘31 

1.050  x  1049 

2.233  x  1019 

2.345 

4.304 

1.535  x  1042 

1.000  x 

1024 

1.535  x  1012 

2.908  x  1041 

2.908  x  10'35 

4.464  x  10'23 

2.233  x  1019 

9.967  x  10' 

3 

1.814  x  1041 

1.000  x 

1024 

1.814  x  1013 

0.000 

0.000 

0.000 

2.233  x  1019 

Total 

2.944  x  1013  Activation  1.679  x  107 
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Appendix  E.  Calculation  of  N20  Steel  Cylinder  Activation 


One  concern  when  irradiating  through  equipment  and  supporting  materials  is  the 
resulting  activation.  The  N20  cylinder  was  to  be  placed  in  direct  line  with  the  beam  port 
so  would  receive  a  full  dose  of  neutrons.  The  activity  induced  could  cause  handling  the 
cylinder  to  be  hazardous  and  require  special  storage  and  disposition  of  the  cylinder. 

Of  particular  interest  was  the  activation  of  58Fe  to  59Fe  with  a  half-life  of  44.5 
days,  which  in  turn  decays  to  59Co.  Iron-58  can  be  gained  from  the  activation  of  57Fe  that 
can  be  gained  from  the  activation  of  56Fe  and  so  forth  all  the  way  back  to  54Fe. 
Calculations  show  that  the  gain  in  Fe  from  the  chain  of  lower  irons  is  negligible  in 
comparison  to  the  amount  of  Fe  initially  on  hand,  even  though  Fe  is  only  0.282%  of 
natural  iron.  An  example  of  this  calculation  done  in  Mathematica  is  shown  at  the  end  of 
this  appendix. 

The  N20  container  was  composed  of  4130  steel  (97.56%  iron)  with  sidewalls  of 
0.5  cm  [9,  32].  The  cylinder  was  20.01  cm  in  diameter  and  51  cm  high.  The  cylinder 
was  placed  upright  against  the  end  of  the  beam  port.  The  beam  was  collimated  to  make  it 
appear  as  a  point  source.  Only  part  of  the  steel  cylinder  was  in  direct  line  of  the  neutron 
beam  as  shown  in  Figure  29. 
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Figure  29.  Geometry  of  cylinder  at  OSU  beam  port. 


The  volume  of  the  cylinder  in  the  direct  path  of  neutrons  can  be  estimated  by 
multiplying  the  thickness  by  the  surface  area  involved  shown  by  the  beam  projections  in 
Figure  29.  The  height  of  the  area  of  interest  at  the  back  of  the  cylinder  is  calculated  from 
using  like  triangles  as  in  Figure  30.  Since  the  height  solved  for  is  half  the  projected 
height,  the  complete  height  is  17.95  cm. 


◄ - ► 


132.43  cm 


Not  to  Scale 


Figure  30.  Using  like  triangles  to  solve  projected  height. 

Vertically  splitting  the  area  of  interest  at  the  shortest  part  and  unrolling  the 
cylinder,  gives  an  estimated  flattened  area  as  shown  in  Figure  3 1 .  This  area  is  then  easily 
calculated  to  give  1043.6  cm"  and  an  associated  volume  of  521.8  cm  . 
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Figure  3 1 .  Estimated  area  of  activation. 

This  area  is  only  an  approximation  of  what  will  be  activated.  Some  scattering 
will  occur  so  areas  outside  the  direct  path  will  receive  some  activation.  The  flux  will  be 
attenuated  somewhat  between  the  front  of  the  cylinder  and  the  back  so  the  activation  will 
not  be  evenly  distributed.  The  approximation  is  only  used  for  getting  an  idea  of  the 
activation  of  cylinder  for  safety  planning. 

Activity  is 


A  =  AN  =  <$>oNT(\-  e~Xt) 


(21) 


where: 


A  =  activity 

.  In  2 

K  =  decay  constant  = - 

halflife 
2  1 

®  =  fluence  rate  [cm'  s'  ] 
a  =  cross  section  [b]  or  [10'  cm'] 

Nt  =  target  nuclei 
t  =  time  of  irradiation  [28]. 

The  thermal  cross-section  (cy  )  is  applied  for  neutron  energies  of  approximately 
0.5  eV  and  lower.  For  higher  energy  values,  the  resonance  integral  (RI)  is  used  to 
represent  the  probability  of  neutron  reactions.  Common  practice  uses  the  RI  just  like  a 
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cross-section  [3].  For  58Fe,  gy  =  1.3  b  and  RI  =  1.2  b  .  The  half-life  of  59Fe  is  44.5  days 

1  O 

or  1068  hours  leading  to  X  =  — — —  =  6.49  x  10'4/h  . 

1068A 

Radiation  will  be  for  7  hours  at  450  kW  (90%  power).  The  number  of  target 
nuclei  is  calculated  from 

N  (22) 

T  AW 

where: 

p  =  density  =  7.874  g  cm'3 

Na  =  Avogadro’s  constant  =  6.0221  x  10“  nuclei  mol' 

AW  =  atomic  weight  =  59  g  mol  1 

V  =  volume  =  52 1 .8  cm3 

%  =  %  of  composition  =  (0.00282)(.9756) . 

This  gives  an  Nj  of  1 . 1 5  x  1 0“~  "  Fe  nuclei.  At  the  end  of  the  205  cm  beam  port 
closest  to  the  core  the  flux  is  7.8  x  10  n  cm'  s'  of  which  4.5  x  10  “  n  cm'  s'  is 
considered  thermal  [1 1].  Of  the  total  flux,  57.7%  is  thermal  and  42.3%  is  fast.  This  is  at 
100%  power.  For  safety  reasons,  the  maximum  operating  level  used  is  90%  so  the  flux 
will  be: 

<bth  =  4.05  x  10 12  n  cm'2  s'1 
(hfast  =  2.97  x  1012  n  cm'2  s'1. 

The  flux  must  be  calculated  at  the  open  end  of  the  beam  port  so,  assuming  no 
scattering  and  a  uni-directional  flux,  the  flux  can  be  divided  by  distance  to  get  an 
estimation.  The  distance  is  205  cm  so: 

<f>th  =  9.64  x  107  n  cm"2  s'1 
Ofast  =  7.07  x  107  n  cm"2  s'1. 
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Then  using  Equation  2 1  and  the  appropriate  cross  section  or  resonance  integral 


gives: 


A*  =  (9.64  x  107  n  cm'2  s''X1.3  x  10'24  cm2)(1.15  x  1023)(1  -  e-(°000649/h><7 
Ath  =  6.53  x  104/s 

Afast  =  (7.07  x  107  n  cm'2  s_1)(1.2  x  10‘24  cm2)(1.15  x  1023)(1  -  e-(°  ooo649/h)(7  ^ 
Afast  =  4.42  x  104/s  . 

This  gives  a  total  activity  of  1.10  x  105/s  or  2.96  mCi. 


Sample  Mathematical  calculations  for  iron  activation  to  59Fe. 


Fe-54  1 

Fe-55  2 

Fe-56  3 

Fe-57  4 

Fe-58  5 

Fe-59  6 

n,  ->  initial  number  of  atoms 


th  ->  number  of  atoms 


cross  section 
d)  ->  flux 


Assumptions:  Flux  is  from  SCALE  at  2xl07.  Decay  from  Fe-55  is  negligible  since  it 
has  a  half-life  of  2.73  years  and  it  is  only  coming  from  Fe-54.  7  hour  radiation  time 
=  25200  seconds.  No  initial  Fe-59. 


DSolve  [ri2  '  [  t] 
{{n2[t]  ->C[1] 


==  ct nil 
1 0  mu.  cti  -i  -i 

1+0  02  i  i 


CT2  +  0*  n2'  [t]  ,  n2[t]  ,  t] 


DSolve [n3 ' [ t] 
{{n3[t]  ->  C  [  1]  + 

DSolve [n4 ' [t] 
{{n4[t]  -+  C [  1]  + 

DSolve [ns ' [ t] 
{{n5[ t]  ->  C[l]  + 

DSolve [n6 ' [t] 


==  ct2  *  0  *  n2  [ t]  -  03*  0  *  (ni3  +  n3 ' [t] )  ,  n3[t] 
r t  -0  ni3  CT3  +  0  ct2  n2  [K$760] 


K$761 

=  =  CT3  * 
t 


dlK$760jj 

1  +  0  CT3  -1  -1 

*  n3  [  t ]  -o4*0*  (ni4  +  n4'[t] 

-0  ni4  ct4  +  0  CT3  n3[K$1417] 


K$1418 


1  +  0  CT4 


n4[t] 
dlK$1417] 


=  =  O' 4  * 
t 


)*n4[t]  -05*0*  (nis  +  ns'ft] 

~0  ni5  as +  0Q4  n4 [K$2096]  dK$20% 

K$2097  1  +  0  05 


t] 


t] 


n5[t],  t] 


==  o5*0*n5[t]  -A*n6[t],  n6[t],  t] 
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{ {n6 [ t ]  ^e"tAC[l]  +  e“tA 


cti  =  2.3*10‘24; 


0  =  2  *  107 ; 
nn  =  2 . 142  *  1024 ; 
ct2  =  13  *  10~24 ; 
a3  =  2.6*  10~24 ; 
ct4  =  2.5*  10~24 ; 
a5  =  1.3*  10‘24; 
ni3  =  3.243*  1025; 
ni4  =  7.358*  1023; 
nis  =  9 . 624  *  1022 ; 
A  =  1.803*  10'7; 


e 

K$2771 


K$2770 A 


0  CT5n5[K$2770]  dK$2770jj 


DSolve [n2 ' [t]  ==  ai *  0  *  nn  -  CT2 *  0  *  n2 ' [t] ,  n2[t],  t] 
{ { n2  [  t  ]  -+  9. 853199999999999  x  107  t  +  C[l]}} 

9 . 853199999999999'*^7  *25200 

2.4830064  x  1012 


DSolve  [n3 ' [t]  ==  02  *  0  *  2 . 4830063999999995 '  *^12  -  ct3  *  0  *  (n^  +  n3  '  [  t] )  , 
n3 [ t ]  ,  t] 

{ { n3 [ t ]  -»  -1. 686359999999355  xl09t  +  C  [  1]  }} 

-1 . 6863599999993546 '  *A9  *  25200  +  3.243*  1025 

3.242999999995751x 1025 


DSolve[n4 ' [t]  ==  ct3  *  0  *  3 . 2429999999957505 ' *^25  -  04  *  0  *  (ni4  +  n4 ' [t] )  , 
n4[t],  t] 

{ {n4[t]  -»  1.64956999999779x10s  t  +  C[l]  }} 

1 . 6495699999977903 '  *A9  *  25200  +  7.358*  1023 

7 . 35800000041569  x  1023 

DSolve  [ns ' [t]  ==  04  *  0  *  7 . 35800000041569 '  *^23  -  05  *  0  *  (n^  +  ns ' [t] )  , 
n5 [ t ]  ,  t] 

{ { n5 [ t ]  ->  3.428776000207846xl07  t  +  C [  1] }} 

3 . 428776000207846 '  *A7  *  25200  +  9.624*  1022 

9. 62400000008 6405 x 1022 

DSolve  [ng ' [t]  ==  05  *  0  *  9 . 624000000086405 '  *^22  -A*ng[t], 
n6[t],  t] 

||n6[t]  ->  1.387820299513292xl013 

2.7  1  8  2  8  1  8  2  8  4  5  9  0  4  51-803xl0"7t  g-1  • 803x10-7 1  +  e-l.803xlCr7  t  c  [i]  jj 
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1 . 3878202995132922 '  *A13  2 . 718281828459045 ' 1 • 8029999999999997  *A~7  *25200 

- 1 . 8029999999999997 '  *A-7  *25200 
e 

1 . 387820299513292  x 1013 

Activity  (Bq): 

1.3878202995132922'*^13  *A 

2.502240000022466xl06 

Activity  (Ci): 

2 . 502240000022466 '  *A6  /  (3 . 7  *  1010) 

0.0000676281081087153 
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